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Abstract 



Quantum entanglement has been recognized as a precious resource in various information pro- 
cessing tasks. Many of its applications rely on the discrepancy between Quantum Mechanics and 
Local Realism. This contrast was first shown by Bell [6J. 

For this reason a large part of the following dissertation has been devoted to the history of 
the discussion on the falsification of Local Realism. We begin with recaUing the consideration 
of Einstein, Podolsky, and Rosen |3l]. Next, we resume the answer of Bohr p2] and the papers 
of Bohm [14l [15], in which he shows a possibility of existence of hidden variables in Quantum 
Mechanics. Finally, we reach the argument of Bell [B], who proved that if such parameters indeed 
exist, they cannot have a local character. 

The Bell argument, originally drawn for two qubits and two measurements per qubit, has 
been generalized to more complicated experimental cases. We present inequalities by Mermin [61] . 
Ardehali [3], Belinskii and Klyshko [5], Werner and Wolf |88], Weinfurter and Zukowski [87], and 
Zukowski and Brukner |103j . All these derivations were done for an arbitrary number of qubits, 
but still with only two alternative measurements on each subsystem. The Bell expressions, which 
utilize more local measurement were presented inter alia by Laskowski, Paterek, Zukowski and 
Brukner [59] . 

Another approach to the problem of finding new Bell inequalities is through an analysis of a 
convex hull, for example, in the space of inter-qubit correlations. This problem has been addressed 
by e.g. Zukowski [102] . His method was further developed for three qubits by Wiesniak, Badzi^g, 
and Zukowski [92]. Therein, we derive versions of the Bell theorem being special forms of the 
criterion presented in [96], but also an inequality, in which all three observers perform one three 
measurements on their particles. Such an inequality cannot be equivalent to any previously known 
ones. 

We also consider consequences of taking the rotational invariance of the correlation function 
as an additional constrain on Local ReaHsm [6^. Upon these assumptions, we introduce a Bell 
inequality, in which every observer performs a continuum of measurements. It is violated by GHZ 
states stronger than any other known inequalities. 

Another problem addressed in this dissertation is the channel transparency necessary to falsify 
Local Realism in the scheme Bjork, Jonsson, and Sanchez-Soto |12| . It is a refinement of the 
proposal of Tan, Walls, and Collett [79], in which a carrier of non-classicality is a single photon. We 
show in [92j that the channel transparency (which in a two-qubit scheme could be associated with 
the quantum efficiency of dtectors) , above which the Clauser-Horne inequality [23| can be violated 
is about 17.2%. This is significantly less than than the threshold on the detection efficiency in 
two-photon experiments, found by Garg and Mermin |37] to be about 82.8% 

We then abandon the problem of the Bell theorem, and pass to the question of entanglement in 
bulk bodies. Because of the high complexity of the physical system and our ignorance about the 
state, quantum correlations must be confirmed with entanglement witnesses [80]. After clarifying 
the notion of a witness we recall some of arguments that low-temperature values of the internal 
energy reveal entanglement [851 fT9] and genuine multi-partite quantum correlations [l^ . We argue 
that also non-linear functions of the state can serve as entanglement witnesses. This is, for example, 
the case of the magnetic susceptibility for systems with rotationally invariant Hamiltonians [90] . 
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The generality of the magnetic susceptibiHty as an entanglement witness is demonstrated with 
various examples. We show that also the heat capacity can reveal entanglement [H]. The last 
result is related to the Third Law of Thermodynamics. 
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Chapter 1 

Introduction 



1.1 Historical Prelude and Motivation 

The formulation of Quantum Mechanics started with introducing the concept of a portion of 
energy in the black body radiation by Planck [60]. Subsequently, Einstein deduced the existence of 
the quantum of Hght [S^ from the behavior of entropy of the electromagnetic field in a cavity when 
one varies its volume. The emerging new theory of Nature, based on wave-functions as a description 
of the system, allowed to explain many phenomena. A lot of attention in the development of the 
theory has been attracted by a phenomenon first discussed in more details by Einstein, Podolsky, 
and Rosen [34j . 

Schrodinger |73] named this unique feature of quantum-mechanical systems entanglement and 
captured its nature (at least for pure states) in an observation that while we may posses the 
maximal knowledge about a composite quantum system treated as a whole, in the extreme case 
measurable quantities related to individual parts of the system may be completely undetermined. 
The subsystems are describes only in a reference to each other. 

Quantum Mechanics allows to arbitrarily superpose wave-functions, even if there exists a 
Nature-preferred basis of distinguishable situations. For instance, an atom can be driven to an 
arbitrary superposition of the ground and excited states. The feature of entanglement directly fol- 
lows from the principle of the superposition. If a joint pure state of two or more quantum systems 
cannot be expressed as a tensor product of the states of each subsystem, but it is necessary to use 
a superposition of such products, such a state is called entangleco ; 

entangled ^ ^ \'4,)W\^)m_ (u) 

A mixed state is entangled if it cannot be decomposed into a convex combination of mixed product 
states, 

entangled ^p^J^ P^pf^ ® ' ■ (1-2) 

i 

i 

The right-hand side of the inequality in l|1.2p is called a separable state. 

Einstein, Podolsky, and Rosen [34] began a discussion on the the possibility of excluding any 
local and reaHstic description of the Universe and posed a question about the completeness of 
Quantum Mechanics. The violation of Local Realism by entangled states was some time later 
elegantly demonstrated by Bell [6j in a form of an inequality. This inequality was thereafter exper- 
imentally confirmed, first by Aspect et al. [4], and generalized for systems of higher complexity, for 
instance, in [6T | [3j \E\ 188 } [87 } [103^ f96 t 159) . or more recently in [102) . This intensive search for new 

^Original german: Verschrankung. 

^Throughout the thesis, the corresponding subsystem shall be denoted by an upper square-bracketed index behind 
mathematical objects, e.g. operators or states. 
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entanglement tests has been motivated by identifying correlations of purely quantum nature as a 
resource for a teleportation of an unknown quantum state [9], quantum computation [42 ^ \77 \ [26]. 
as well as decreasing the communication complexity in computational tasks [21]. It has been also 
assigned a role in certain protocols for distributing a secret cryptographic key [51 [Ml ES] (The 
relation between entanglement and communication complexity problems or qunatum key distribu- 
tion is through the violation of Bell inequalities). However, protocols for the last two tasks were 
proposed and experimentally realized also with a use of a single photon [32l [H3 ■ 

However, even if an entangled state does not violate Local Realism, it can be useful from the 
point of of view of quantum information processing. It was shown by Bennett et al. [10] that 
weakly entangled states can be locally post-processed with local actions classical communication 
communication between partners possessing the state (in many copies) in order to distill maximally 
entangled pairs of qubits. Even if it is not possible for a given entangled state, that is the state 
is bound entangled [49], its correlations can be used to activate non-classicality of other physical 
systems [17]. In any case, non-separability of a state can be confirmed with positive, but not com- 
pletely positive maps [67\ I49j or entanglement witnesses [801 |54]-(originally) hermitian operators, 
extreme eigenvalues of which are associated to entangled states. 

1.2 Summary of Results 

For the reasons explained in the first section, a large part of the research, in which the Author 
of the following dissertation has participated over last years, has been devoted to the problem of 
new versions of the Bell Theorem. Results with potentially practical applications are presented in 
Chapter 3. Therein we exemplify a method of generating new Bell expression from the analysis 
of a convex hull in the correlation space for the simplest non-trivial case of three qubit^ and 
a choice of up to three observables for each observer. For most of the newly-found inequalities, 
we present a sufficient condition on a state for satisfying the inequality. The Chapter ends with 
considering consequences of taking the rotational invariance as an additional constraint on Local 
Realism. Assuming this we are able to present a Bell condition, that can be violated exponentially 
stronger by GHZ states than any thus far known standard inequalities. 

In Chapter 4 we argue that the so-called "single-photon" Bell-type experiment, initially proposed 
by Tan, Walls, and Collett [T^ and theoretically refined by Bjork, Jonsonn, and Sanchez-Soto ^12j . 
is very robust against photon losses. 

Entanglement witnesses allow to investigate non-separability in bulk solid systems, arbitrarily 
large lattices of spins in a thermal equilibrium. We then aim to detect quantum correlations 
with certain values of thermodynamical quantities, like the internal energy, the heat capacity, 
and the magnetic susceptibility. The pioneering work due to Wang and Zanardi [85] Hnks the 
internal energy with a measure of two-qubit entanglement, the Wootters concurrence [94J- In spite 
of possible experimental difficulties in measuring the internal energy, many authors [HI [Ml [ISl 
(951 [44] followed the original idea. On the other hand, Brukner, Vedral, and Zeilinger [19] have 
shown that entanglement is implied by low-temperature values of the magnetic susceptibility for a 
specific substance, cupric nitrate (Cu(N03)2)- In section 5.6 we will repeat the argument from [90] 
that the magnetic susceptibility is a valid entanglement criterion for thus far the widest class of 
materials, all those in which spin interaction is isotropic. Also in Section 5.8 we aim to prove that 
independently from symmetries possessed by the Hamiltonian, entanglement can be revealed by 
low-temperature behavior of the heat capacity ^1] . We stress that this result is strongly related to 
the unattainability of the the absolute zero temperature, which is possibly equivalent to the Third 
Law of Thermodynamics |65| . 



qubit (quantum bit) is a physical system described by a two-dimensional Hilbert space, for example a spin-^ 
or a polarization of a photon, as well as a unit amount of quantum information, that can be stored in such a system. 
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Chapter 2 

Bell Theorem-Introduction 



2.1 Einstein-Podolsky- Rosen Paradox 



In 1935 Einstein, Podolsky, and Rosen |34j noticed the discrepancy between the point of view 
that Quantum Mechanics is a compete theory and the probable from their point of view description 
of the Universe based on two main assumptions: 

• Locality is defined from Special Relativity [33] • It states that no two events can infiuence 
each other if they are spacelike separated. The mutual independence not necessarily means 
that the events might not have had the common cause. 

• The other important notion, a new meaning of which was introduced in [34], is Realism. The 
Authors themselves, avoiding complicated philosophical considerations, conclude that: 

If, without in any way disturbing a system, we can predict with certainty (i.e., 
with probability equal to unity) the value of a physical quantity, then there exists 
an element of physical reality corresponding to this physical quantity. 

Thus Realism means that the results do not appear at the moment of the measurement, but 
rather wait to be revealed by an observer. Such a requirement is natural in the macroscopic 
world. For example, objects have their dimensions and masses, although they could have 
never been measured or weighted. The act of the measurement only reveals facts unknown 
to the observer. 

The theories based on these two assumptions are called local reaHstic theories or theories with 
local hidden variables (LHV). Local Realism can be interpreted as a situation in which every part 
of a quantum system carries a set of information about which result would be yielded under any 
measurement. 

Non-avoidable randomness in Quantum Mechanic forced the Authors to face the following 
dilemma: either (i) the quantum-mechanical description of the Universe based on a wave-function 
cannot be considered complete, that is not all the elements of the reality are described, or (ii) 
quantities related to non-commuting observables cannot be simultaneous elements of the same 
reality. 

The potential paradox was illustrated with the following wave-function of two particles with 



where xq is a constant and p here denotes a momentum of either of particles. The integration gives 



coordinates x'^' and x^^^ 




(2.1) 



^(xW,2:''')oc(5(a;W-a;l2l-Ka;o), 



(2.2) 



^It is enough to consider particles living in one-dimensional geometric space. 
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and after the double Fourier Transform we get 

^ + p[2] _ (2.3) 

with po being another constant. The particles are thus separated by a constant distance, and in 
the same time the sum of their momenta is equal to po- As noticed by Schrodinger [73], when we 
trace out one of the particles it turns out that the other particle is completely delocalized, in both 
the position and the momentum representations. Particles can be equally likely found in every 
point of the Universe and their mean kinetic energy is infinite. 

According to l|2.2p if one of the observers (hereafter called Alice) detects a particle at some 
given point x, she and only she knows instantaneously that the second particle can be found at 
xq + X. At the same time the other observer (called Bob) would find it exactly there in spite of 
the fact that due to the assumption of locality, no information about the outcome of the Alice's 
measurement could have propagated with superluminal speed. Thus, following the definition, the 
position of Bob's particle is an element of the reality. But Alice could have decided to measure the 
momentum of her particle rather than position and get some value p and learn that Bob's particle 
is has the momentum —p + po, which also seems to be an element of the reality. 

One must remember, however, that operator x ^^^^^1 does not commute with p I^^^)]^ g^j^j ^jj^g 
Quantum Mechanics forbids to have well-defined values of the position and the momentum of the 
same object at the same time. If we agree that both the position and the momentum of the second 
particle are elements of the same reality, the Quantum Theory cannot be complete, as according to 
this theory only one of the values can be fully defined at a time. If we insist that a wave-function is 
the most complete description of a system, the position and the momentum cannot be simultaneous 
elements of the reality. On the other hand, different representations of the same wave-function are 
related to each other by transforms. Einstein, Podolsky, and Rosen had come to the conclusion 
that Quantum Mechanics could not have been considered complete. 

2.2 Anwser of Bohr 

A direct response to the EPR letter pi| was an article of Bohr [16]. Using the notion of 
complementarity, which was introduced by him, he explained that considered criterion of com- 
pleteness of Quantum Mechanics is not adequate due to fundamental differences between classical 
and quantum-mechanical descriptions of the reality. 

The paper starts with introducing two theoretical situations. In the first one, a particle falls 
on a diaphragm with a slit and, provided that it was aimed at the slit, its trajectory is distracted 
by a diffraction. The diaphragm is a part of a larger measuring device and is fixed with respect to 
the coordinate system of the laboratory. The only uncertainty of the position of the particle just 
behind the slit is related to its width, here denoted as Aq{— ^ {cp-) — (g)^, (•) denotes the mean 
value). The diffraction causes a momentum interchange between the particle and the diaphragm 
of the magnitude Ap(= V (p^) — (p)^)- This momentum is, however, transported to the optical 
table, and further to the laboratory. It is hence impossible to precisely determine Ap. Therefore 
we cannot predict the exact point of a far screen at which the photon would be detected. 

In contrast, we can consider the case in which the diaphragm is a freely movable part. Before 
and after the photon passes the slit, we can measure the momentum of the diaphragm with help 
of a stream of probe particles. Again, any such measurement would mean a displacement of 
the diaphragm due to coUisions with probe beams. Although we have observed the momentum 
interchange, we deal with the uncertainty of the position. Still, we are unable to predict the exact 
point of the detection. 

How does this correspond to the EPR pair? Instead of one, let us consider two slits in a 
diaphragm which are both narrow in comparison to the distance between them. In case of movable 
slits, when we are able to measure momentum of the diaphragm, the state of the two particles 
is close to the EPR state. We can obtain knowledge of the difference of positions jg^^l — as 
precisely as narrow are the slits, and, under certain assumptions, we can learn about the amount 
of the momentum interchanged with an arbitrary precision. However, since the initial position of 
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the plate was unknown, we cannot determine g'^l + gl^' , and as it was a three-body collision, it is 
also impossible to have a precise value of Now it is up to the observer which quantity he 

measures (on one particle) and which corresponding property of the second particle will be defined. 

The difference, from the point of view of Bohr, lies in the context of the measurement. In the 
first case, the diaphragm was a fixed element of a measuring device. We can arbitrarily precisely 
determine the position of the particle in the initial state. In case of a movable diaphragm, it 
was rather a part of a quantum system. In a sense, a single slit experiment is equivalent to the 
EPR gedanken-experiment, with the diaphragm playing a role of the other particle. In both cases 
Quantum Mechanics allows to predict certain properties of the particle without any disturbance, 
thus it cannot be considered incomplete. The division of the system to a microscopical quantum 
system and a macroscopic measuring device is the basis of Bohr's Copenhagen interpretation. 
Significantly, we are always forced to deal with an unavoidable and uncontrollable interchange of 
certain quantities, like the momentum, between the two systems. 



2.3 Bohmian Interpretation of Quantum Mechanics 

The discussion on completeness of Quantum Mechanics returned on the occasion of papers by 
Bohm |141ll5j from 1952, in which he proposed his own interpretation of the theory. He argued that 
the formulation of the Quantum Theory is fully consistent and it is the most possible complete 
theory of nature, however it is still possible to introduce additional parameters. They are, by 
definition, not accessible for the observer. These parameters shall allow to predict results of all 
possible measurements, as Einstein, Podolsky, and Rosen were considering. He based his statement 
on the analogy that Thermodynamics correctly describes the behavior of gases, although it neglects 
positions and momenta of individual atoms and just gives their statistical distributions. Similarly, 
Quantum Mechanics may contain some elements, which carry information about outcomes of all 
possible measurements, but since they are never known, the theory gains a statistical nature. 

The Bohmian interpretation is based on the new understanding of the Schrodinger equation. If 
we write the wave-function ip{x) of a particle of mass m in a potential V{x) using two real-valued 
functions ipix) = i?exp(iS'/ft), the Schrodinger equation, 

= -|^VV(x) + V{xMx), (2.4) 

no 1 

ds {{ysf , v^R\ , , 



can be expressed as 



dt V 2m ^ ' 2m R 

It is convenient to introduce the probability density P(x) — R^{x): 

Then, in the classical limit of ft — + 0, the phase function S becomes a solution of the Hamilton- 
Jacobi equation, so that we can interpret its gradient as the momentum. I|2.7p expresses the 
probability conservation. We thus see that that the Schrodinger evolution differs from the classical 
one in such a way that besides the "classical" potential we have the "quantum-mechanical" one, 

^ 4m V 2P^ J 2m R ^ ' 

This greatly resembles the situation in Electrodynamics. In both cases the field, the electro- 
magnetic field or the wave-function satisfy certain equations (Maxwell's or Schrodinger's), from 
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which we compute forces. Thus knowing the state at a certain instant of time we are able to trace 
the whole evolution, back forth and back in time. 

The essential difference between the electromagnetic field and the wave-function is the form of 
the equation. The Schrodinger equation is homogeneous with respect to "0, whereas the Maxwell 
equations are not with respect to the field. Bohm sees this difference as a possibility to modify 
Quantum Mechanics, such that significant effects would appear only at very short distances. 

The inhomogeneous modification of the Schrodinger equation suggested Bohm emerged from 
the fact that at the time some phenomena in the subatomic scale did not have proper descriptions. 
Quantum Mechanics might turn out to not be valid at the Planck scale, that is distances of the order 
of 10~^^ meters and would then need to be modified. Obviously, a more general theory, correct also 
at these distances, must reproduce results of Quantum Mechanics. Bohm now suggests to abandon 
his idea of seeing VS'/m as a velocity v of the particle, as well as to assume that the homogeneous 
and linear with respect to ip Schrodinger equation is not an evolution equation. Giving up these 
two hypotheses implies that we cannot treat the wave-function as a statistical description of an 
ensemble of particles. 

An example of such a correction in the papers of Bohm is done by introducing a force, which 
would equahze the physical velocity w of a particle and VS'/m. This should happen in a very short 
time T of the order or 10~^^m/(|l or 10~^'^s. A modified classical equation of motion has a form of 

and the non-homogeneous Schrodinger equation reads 

ih^^H^/j + ap-^S), (2.11) 

where /(O) = ^(0) = 0. The quantum evolution returns to its homogeneous form after time r. 
Until then, all basic principles of Quantum Mechanics, such as the superposition principle, might 
not be valid. However, this or other corrections of such a character would allow to extend the 
theory with hidden variables predicting simultaneously the position and the momentum with an 
arbitrarily high precision. 

In [15] the Author explains his interpretation of the Einstein-Podolsky-Rosen gedankenexper- 
iment. According to Bohm, since the wave-function l|2.2p is real, both particles remain at rest. 
Their positions are described by a probability distribution, such that always x'^^^ — x'^'^^ = xq. The 
observer measures the position of one of the particles and disturbs the momentum of the other in 
an uncontrollable way, and, similarly, the act of measuring the momentum of one particle unpre- 
dictably disturbs the position of the other. In this way Bohm is able to save the statement that 
one cannot measure two arbitrary quantities, even though in case of the EPR experiment one acts 
on two spatially separated subsystems. The disturbance of field S may happen with a superlumi- 
nal speed. This hypothesis must be of course questioned with respect to its agreement with the 
Special Relativity, as predictions of the latter are not in a confiict with Quantum Mechanics. An 
instant change of the quantum- mechanical potential cannot result in transfer of any useful informa- 
tion between two distant observers. Without a measurement neither of them has any information 
about the position or the momentum of his particle until he gets a classical information about the 
other half of the EPR pair. The transfer of the classical information only via the EPR channel 
would be possible if there was a way to predict the position or the momentum without performing 
any measurement. This is clearly impossible in Schrodinger's formulation of Quantum Mechanics, 
which, as Bohm argues, does not hold at Planck scales. One can construct theories, which, as it 
was shown, allow to overcome the Heisenberg bound for uncertainties, but also these, which permit 
superluminal propagation of informations. Then, we have two more options to avoid a discrepancy 
between Quantum Mechanics and Special Relativity. Firstly, there might exist a law, which simply 
prohibits superluminal changes of the controllable parts of the quantum-mechanical potential also 
below the Planck length. The other option is that the Lorenz invariance can be just irrelevant at 
nanoscales. It would be an important extension of the ideas of General Relativity to consider the 

here denotes meters, not the mass. 
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metric tensor g, and thus the curvature of the spacetime, dependent on the wave-function ip{x), 
which describes the matter and the energy in this space. The postulate of the speed of Hght as the 
universal constant, the largest possible speed in the Universe, could then be relaxed. 



2.4 Bell Theorem 



Bohm had argued [14tll5). that Quantum Mechanics under some assumptions can be interpreted 
as a deterministic theory based on hidden variables. However the paper of Bell [6| from 1964 showed 
that these additional parameters, if they exist, cannot have a local character. 

Bell proposes that AHce and Bob can locally perform one of two possible dichotomic measure- 
ments, A[^^,^2^' for Alice, A!j^^,A2^' for Bob, which yield "-I-1" or " — 1" as results. Let us denote 
the results as af', 03^' for AHce and a^^', Oj^' for Bob. The assumption of Realism allows us to use 
all four quantities in the same time and Locality excludes mutual influences between events at 

Alice's and Bob's side. Out of two expressions, a!j^^ -I- Oj^' and af' — Oj^', exactly one is equal to 

[21 [21 [21 [21 

±2, the other vanishes and similarly for + and — ■ As a consequence, for any given 
set of leal hidden variables determining results {c^^\ a\^\ ol\^\ a^2^} , only one of four expressions 
(aW±aW)(afl±4'l) is equE 
the experiment we can write 



(a^^' ± a2^')(a5^' ± Oj^') is equal to ±4, while three others are 0. After averaging over many runs of 



s[ll,s[21=±l 



(aW + ,s[iUW)(421 + s[2U|,21; 



<4, (2.12) 



where the subscript A stresses that l|2.12p is vaHd for a given set A of local hidden variables. 

One may also consider theories, in which various possible sets of LHV are distributed with a 
probability density p(A). This is reflected in our uncertainty about measurement outcomes. To 
compute averages in such stochastic theories one must integrate over all possible sets A with p{\) 
as an integration kernel: 



E 

s[il,sPl=±l 



(aW+,sWaW)(4^1+sP14^1)\ p(\)d\ 



X 



< 4, (2.13) 



where dX is a metric in the set of A. 

Since the Bell inequality is linear with respect to the mean values, its local realistic bound 
can be concluded from the flrst kind of theories. We hereafter assume that the theory under our 
consideration is deterministic and denote local realistic predictions by a subscript LHV. 

Since taking the modulus changes at most the global sign of the expression, l|2.12p can be also 
written with use of a sign function S{s^^^ = ±1, s'^' = ±1) = ±1 [103J: 



^(S[1U[2])^(^[1]^^[1]^[1])(^[2]^^[2]^[2])\^^^ (2.14) 



4< E 

s[il,sPl=±l 



From the point of view of the discussion on the possibility of local hidden variables only these 
eight inequalities l|2.14p are important, in which the sign function is not factorisable, S'(s[^l, s'^') ^ 

S'i(sI'^l)S'2(s[^l). The remaining eight produce trivial bounds (y4[^'4^') ^ 1- In the flrst case, the 

inequalities are equivalent to the one introduced by Clauser, Horne, Shimony, and Holt [24|: 

- 2 < (aW^I^I + 4114^1 + 4Ufl - 4'4">LW < 2. (2.15) 

These inequalities may not be valid anymore when one uses the formalism of Quantum Me- 
chanics to compute mean values of measurements. Local measurements with results "±1" on 
two-dimensional quantum systems are associated with observables in a form of a scalar product of 
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+1 -1 




+1 -1 




^Qur^ 










Alice 




Bob 



Figure 2.1: A scheme of a Bell experiment. Two partners, Alice and Bob, receive particles from 
a common source and perform one of two measurements of their choice. The mean values are 
plugged into ((2A21 1 



a normalized three-dimensional vector built of real components dj with a vector of Pauli matrices 

"'^{l J ) ' = ( ° V ) ' ^3 ^ ( -1 ) ' (2.17) 

so that Alice's observables are given by Af^ = di • a and Bob's-^f — • (jl^l . The the 
eigenstates of the operator are represented by |0) ^ ^ ^ ^ and |1) ^ ^ ^ ^ , respectively. 

Now, let us assume that the pair of spins— ^ is in the singlet state 1^*^) = (|10) — |01))/\/2, 
where in the state |0) the spin is parallel to the z-axis and anti-parallel in |1). Then the correlation 
function is given by (Af^A^^^) = — -a^^' . For convenience, we will restrict ourselves to observables 
in the xy-plane, which are given by 

= cos(^ACTr' +sin(p^aW = fvA o''^ ) ( 1 ) ' ^^'^^^ 

A^'^i^B) = cos^saf^ + sin(pB4" = ( J l) ^ { ' ) " (^'^^^ 

The correlation function is then equal to 

E{^^a.^b) = (*-|AW((/?A)APl((^i3)|*-) = -cos(^a-<^b). (2.20) 

Let us now choose two observables at each side. It turns out that the optimal choices are = 
AW(0),4^1 = AW(7r/2), ^[^1 = A[2l(7r/4),A|1 = A[2l(-7r/4). The middle of (l215ll takes then the 
value of — 3cos7r/4 + cos37r/4 = —2-^2. This is less than any theory with local hidden variables 
can predict. The Bell theorem thus states that there exist measurements, for which quantum 
mechanical predictions are in disagreement with all possible predictions based on Local ReaHsm. 
The natural question arises whether this can be confirmed experimentally. 



2.5 Additional Assumptions for Bell Theorem 



^These observables, together with (tq = ^ ^ ^ will appear often in our further considerations. They con- 
stitute an orthogonal basis of 2 x 2 matrices, as Traiaj = 25ij. It turns out that it is extremely convenient to 
parametrize an Af-qubit state by mean values of products of Pauli matrices, 

Ta...n = Trp<7W...<7|fi; p = ^ (El....„=o7^-.."'^™-'^n'^') ■ (2-16) 

The set of numbers {Ta...n}a „=o shall be called a correlation tensor. Strictly speaking, only {Ta...n}a n=i 
transforms properly under 0(3) rotations, and originally the name has been given to this subset. The rest of 
coefficients describes reduced states. This two meanings of the notion of the correlation tensor can be met in this 
dissertation. 
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Figure 2.2: The experiment of Weihs et al. [86|. The intervals between generations of random 
numbers by RND is about 10 times shorter than the time of the propagation of photons in fibers. 

Apart from the Locality and Realism we need to take one more very important, but often 
forgotten assumption. We need to assume that the observers choose local observables not long 
before the measurements and independently from each other. This could be guaranteed if both 
of them use fast random (or even pseudo-random) number generators to determine which local 
measurement shall be performed. If the choice of measurements for each run of the experiment 
was predefined much before the emission of the entangled pair, a local reaHst could claim that the 
hidden parameters of the particles were infiuenced by this predetermination. 

The conditions of Randomness and Locality were strictly satisfied in the experiment performed 
by Weihs et al. [86]. In this realization of the experiment, each measuring station was connected 
with a source by 500 meters of an optical fiber, but one half of each cable was rolled up next to the 
source, while the rest was strechted to AHce and Bob, separated by about 400 meters. Entangled 
pairs remained for about 1 ^s close to the source, and propagated to the observers for roughly the 
same time. During the propagation time local measurements were randomly chosen using ultra-fast 
opto-electronic switches. Random number generators consisted of a light emitting diode, a filter, 
a beam splitter, and two detectors behind it. The filters attenuated the emitted light to about one 
photon per pulse, which was detected at one of the sides behind the beam splitter. The results of 
these measurements determined the position of an opto-electronic transducer, through which an 
entangled photon passed, thus the choice of a local observable. Such a decision was made 10'' times 
per second, finally when photons were not further than 25 m from the stations. Information which 
observables were chosen could have reached in time neither the source, nor the second particle. 

In spite of rigorously satisfying the Locality and Randomness conditions, the experiment of 
Weihs et al. [8^ was not ultimately convincing, since the quantum detection efficiency, that is 
the probability that the particle was detected given that it actually reached the detector, was 
insufficiently high. 

Let us denote the detection efficiency of every detector used in the experiment as rj. Following 
Garg and Mermin [37] , for the sake of view of advocating Local Realism it is convenient to assume 
that each non-registered photon would produce always the same result, (say, "-t-1"). Let us, 
moreover, assume that the state is a Werner state, that is a mixture of a maximally entangled 
state, |«'~)(*"| = 5(|01) - |10))((01| - (10|), with a maximally mixed state, ^^fl: 

p = V\^-){^-\ + {l-V)'^, (2.21) 
< Z < 1. 

Such states are created when one of the photons from an entangled pair is being a subject to a 
depolarizing channel described by a transformation p Vp + {1 — V) J2'i=o ^l^'pf!^' ■ Thus after 
taking into account the detection efficiency, we have a effective correlation function E^ffi^A, ^b) — 
—Vrf cos^LpA — 'Pb) + (1 + ■ After plugging into l|2.15p . one gets 

- 2 < -VTf{co&{pA - Pb) + cos((^^ - pb) + 

cos(^A - p'b) - cos(^^ - p'b)) + 2(1 - 7]f < 2, (2.22) 

■^By ixxx we denote a unit matrix of dimension x. 
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and optimizing over angles we obtain a necessary condition for a violation of the inequality: 



Vri^V2 + {l-r]f >l. (2.23) 

This is a condition for critical values of V and rj only above which the Bell inequality can be 
violated: 

^•^^^^ = /OV ^ A-V (^-^^^ 
\ 'iVCRIT + i 

The inequality (|2.15p cannot be violated if ^ < ■ If the photons reach detectors without any loss 

of coherence, the required detection efHciency to violate a CHSH inequality is tjcrit — 2(^/2— 1) « 
82.8%. 

The relation (|2.24p can be also obtained from the Clauser-Horne (CH) inequality [25 for event 
probabilities: 

PiipA, iPB) + P{v'a. Vb) + P{^A, v'b) - P{^'a, v'b) - PiVA) - P{vb) < 0. 

(2.25) 

2 

The joint probabilities take a form P(iy9A, V-b) = ^{^~^ '^'^K'Pa — ^b)), whereas local probabilities 
are given by P{ipa) = P{^b) — §• ^PPlying this forms to l|2.25p and after the optimization we 
get ^ 

!L(i_V21/)-ry<0, (2.26) 

and finally obtain (|2.24p . 

It was shown by Eberhard [3T] that the critical efficiency can be lowered to 66.7%. This is 
done, however, for states arbitrarily close to product ones. 

One should also consider false detections, called dark counts. As we assume, they occur with 
the probability Pd within the time gate. Now, the effective probabilities are Peff{(pA,y^B) = 

^(1 - VcOs{ipA - ipB)) + »7(1 - "i)PD + (1 - §)'P|, and Peffi^A) = Pi^PB) = i + (1 - ^)Pd, 

which put into l|2.25p gives 

VcRiT = (1 - Pd){2 - T]cRiT){PDi2 - t^crit) + TlCRiT)l{V2rfcRiT)- (2.27) 

The dark counts have their origin in thermal excitations of electrons in a photosensitive element, 
which are not distinguishable from photoexcitations. To avoid this problematic effect, detectors 
are often cooled down to the temperature of liquid nitrogen, i.e., 72K. As a controllable parameter, 
the dark count probability is not as relevant as the detection efficiency. 

Out of many realizations of the Bell experiment, only in the one of Rowe et al. [71] the detection 
loophole was overcome. They have realized an entangled state with trapped ions ^Be"^, what has 
allowed to achieve almost perfect measurement efHciency (about 99%). However, the two entangled 
ions were placed in the same trap, and separated from each other only by micrometers. Thus the 
locality loophole was left wide open. 
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Chapter 3 

Generalizations of Bell Theorem 



3.1 GHZ paradox 

In 1989 the two-qubit Bell theorem [6] had already been broadly discussed in literature (e. g. 
[231 im [8]), and there have been realizations a Bell-type experiment [1] by then. Although up to 
date no such experiment had all the possible loopholes closed, already at the time the majority 
of the physical community accepted the possibility of having two quantum objects correlated in 
a strictly non-classical way. One natural question was if entanglement can be demonstrated with 
more than two mutually correlated objects. Another interesting problem was to concentrate the 
research on systems described by Hilbert spaces of dimensions higher than two. At the first sight, 
both Hmits, the Hmit of large numbers of subsystems and the one of many degrees of freedom, may 
suggest reaching macroscopic, thus fully classical regime. However, it is just enough to consider 
two spin-^ particles, which isotropically interact with each other, 

i? = 5W-5[2l. (3.1) 

The ground state is a singlet state, in which no information about individual spins is available. 
The von Neumann entropy of the reduced states, which since 1996 is known to be a unique measure 
of entanglement for pure, bipartite systems [7], reads 

5(Tr2|ff)(ff|) - -Tr(Tr2|g)(5|)log2(Tr2|.g)(5|) 
= -Tt^^(H!±i)^ log, ^ ^ ^) (3_2) 

and grows logarithmically with the magnitude of the spin. Here Tr2 denotes a partial trace over 
degrees of freedom of the second particle. One can also show, for example using an argument from 
[ST] or [90], that the higher the spins are, the higher is the critical temperature, above which known 
criteria do not reveal entanglement. Can more qubits also lead to more non-classical predictions 
than two qubits? 

The first positive answer to this question was given by Greenberger, Horne, and Zeilinger 
|41) . who have exemplified an even stronger discrepancy between predictions of Local Realism 
and Quantum Mechanics than found in [6j. Surprisingly, this did not require a new inequality. 
Their original derivation was done for four qubits, nevertheless, let us discuss the version for three 
particles in more details. 

Let each of three observers, Alice, Bob, and Charlie, receive a spin-^ coming from a common 

source. Each of them chooses one of two dichotomic observables (A ^ , ^ ^ , ylf ' , vl'^l ,Af\ A^^^ ) 
with spectra {-1-1, — !}• Let these observables be projections of the spins onto axes in the xy-plane, 
as in the case of two-qubit Bell theorem. Assume that the results of these measurements were 

^The absolute zero temperature is not attainable in a finite process |65| . but temperature can be made arbitrarily 
close and hence the actual thermal state can arbitrarily well approximate the ground state 
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deterministically predicted by local hidden variables. For example local hidden variables may set 
three different products, U'^U^^' , A^^Uf^U^^' , A^^U^^U^^' to be equal 1. Then from these three 
results we can compute the result for Aij^'^J^' : 



Af Af = 1 

Af Af Af = 1 

Af Af Af = 1 

Af Af Af = 1 



(3.3) 



This is due to the fact, that since the spectra of all observables are {+1,-1}, their squares are 
just unity operators. On the other hand, if AHce, Bob, and Charlie measure x or y components of 
their spins, for a certain quantum mechanical state they can have 

iafaf.f) =1, 
iafafaf) =1, 

iafaf.f) =1, 

i'rfaf.f) =-1, 



(3.4) 



what follows from 



af af af 
af af af 
af af af 



X Jl] _^[2]^[3] ■ (3.5) 



a, —a, CTi 



A quantum state, for which mean values are equal to those assumed in l|3.4p takes a form 
of -^(lOOO) — I 111)) and is referred to as the 3-qubit GHZ state. Similar calculations hold for 
other possible results of first three products. The calculations based on Local ReaHsm give results 
opposite to quantum-mechanical ones. 

This way of showing the contrast between Local Realism and Quantum Mechanics is called "All- 
versus-Nothing" Bell theorem, or the Bell theorem without inequalities. It can be easily generalized 
to more qubits. A possible set of products, which allow to make predictions for all non-vanishing 
mean values in xy-planes consists of HiLi <^i^ and all the combinations, in which AHce and one other 

observer choose a2, while the rest of them measures ci, e.g. afaf JliLs '''i'' 7 '''2^' '''2^' 11^=2 4 '^1^' 
etc.. 

All-versus-Nothing paradoxes can be also derived for hyperentangled H pairs of photons [22]. 

It is also possible to derive All-versus-Nothing paradoxes for higher dimensional systems. For 
example, for three qutrit^ Alice, Bob, and CharHe choose one of two observables with spectra 
{1,Q!, a^} (a = exp(27ri/3)), so that (^P')'^ — t^xaii — l,2,j — 1,2,3). In such a case, a local 
realistic product 



Af 


Af 


Af 


Af 


Af 


Af 


Af 


Af 


Af 


Af 


Af 


Af 


Af 


Af 


Af 



(3.6) 



^By hyperentangled we mean entangled in more a pair of degrees of freedom, such photons are entangled both in 
polarization modes and spatial modes referring to different paths in Mach-Zehnder interferometers, in which they 
propagate. 

^Systems of the Hilbert space dimensionality 3. 
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is relevant. If we find a set of observables and a state for which, e. g., = 
(4''^f'^?'> = {A^l^AfAf) = (aWaI'U'^I) = 1, but (aW^'U^'I) 7^ 1, one has an explicit 
violation of Local Realism. Nevertheless, paradoxes derived for higher-dimensional systems are 
usually not as strong as those for qubits. 



3.2 Mermin Inequalities 



Mermin (6^ observed that the whole non-classicality of the 7V-qubit state 

|GiJZAr) = -i=(|0)«^ + |l)«^), (3.7) 



the correlation function in xy-planes of which reads 

N 



E{(^^,...,^n) = (gHZn 



J^(cos0i(T[ + sin^icr. 

i=l 

is related to the fact that it is an eigenstate of the operator 



GHzJ^= cos (Y.^^^ ' 



(3.8) 



/ N N 



M=-{ ||«+*CT^^')+||(4^1-^4')|- (3.9) 



In l|3.9p only these terms survive in which (72 appears /c = 0,2,4, ..,2/ni(7V/2j3times. If/c/2 IS even, 
there is a sign in front of a product of Pauli matrices, and " — " otherwise. The corresponding 
mean values for the GHZ state (|3.7p have the same sings and each of them has a modulo 1. Thus 
the eigenvalue of p.9p corresponding to l|3.7p is equal to the number of non- vanishing terms in M, 
2^-1. 

On the other hand, if we assume Local Realism, the mean value of A can be written as 

F - {M)lhv = Re j n(«i'' +*a2 ') ) • (3-10) 

V = ^ J LHV 

By ap' we here denote the outcome of the af measurement determined by LHV. F is a real part 
of a product of N complex numbers, each of the modulo \/2 and the argument ±7r/4 ± tt. The 
vector representing in a complex plane this product can be parallel to the real axis for N even, 
or create angle 7r/4 with the axis for TV odd. The classical limit of this expression is hence 2^/^ 
and 2^^"^^/^, respectively. Thus Mermin |6T| gave the first quantitative argument that the GHZ 
state for N qubits leads not only to exponentially many (with TV) All-versus-Nothing paradoxes, 
but also to an exponentially strong violation of Bell inequalities. 

3.3 Ardehali Inequalities 

The first generalization of the Mermin inequalities was introduced by ArdehaH \3} . He suggested 



4 



By Int{x) let us here denote the integer part of x, e.g. Int{2.5) = 2. 
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to consider two operators for the A^-qubit case: 



Ml 




+ ...) 



+ ...) 



(3.11) 



M2 




[AT-l] 



+ 



(3.12) 



As before, the scalar product of real unit vectors a [ and a 2 with the Pauli matrices vector 
defines the first and the second observable of the TVth observer. The first operator is built only of 
terms in which 1T2 appears k — 0,2,4, ...,2Int{{N — l)/2) times for qubits 1 to iV— 1. Is fc/2 even, 
the sign in front of the term is " + " rather than " — ", which appears for fc/2 odd. Similarly, in 
the second operator (T2 appears fc = 1,3, ...,2Int{{N — l)/2) + 1 times within this set of qubits. 
There is a " + " sign for (fc — l)/2 even and we have " — " for (fc — l)/2 odd. Thus the operators 
can be written as 



Quantum mechanically, if vectors a and 0. 2 chosen to be (1/V2, l/-\/2, 0) and 

(l/\/2,-l/V2,0), respectively, we can reconstruct the original result of Mermin, that is that 
l|3.7p can be made an eigenstate of Mi + M2 with the respective eigenvalue 2^~^/^. 

Local and reaHstic calculations are also similar to the ones made by Mermin. It is enough 
to notice that, as in case of the original derivation, for local realistic theories only one of two 
expressions, (A^^^) LHV + {Ai2^^) lhv and (A^-^^) lhv — {^-^2^^) lhv , is equal to ±2, the other vanishes. 
It then suffices to repeat the bounding argument of Mermin for l|3.13p and p.l4p to obtain that 
the local realistic bound is 2^^"^)/^ for N odd, or 2^^^ for TV even. 

Thus in the case of the GHZ states ArdehaH's inequalities can be violated with a strength up to 
V2 times higher than Mermin's inequalities. This advantage was achieved in [3] by giving the last 
observer a possibility to measure projections of the spin onto two arbitrary axes, rather than onto 
two fixed ones. This allows Adrehali inequalities to be violated by both the iV-partite GHZ state 
and a product of the {N — l)-qubit GHZ state and a single qubit pure state. The derivation of 
the inequality also refers to Bell's original idea of pairs of in some sense complementary quantities. 
The complementarity is expressed in the fact that if one of the expressions takes its maximal value, 
the other vanishes. 

3.4 Mermin- Ardehali-Belinskii-Klyshko Inequalities 

The idea of such a complementarity was exploited in the inequalities derived by BeHnskii and 
Klyshko [5|. Let us consider two series of operators, , S2 , defined in an iterative way: 



To initialize this iteration, we need to assume that the seeds are 

^2^'. The structure of S2 is the same as of , but with A^^^ and ^j'^' interchanged for all k < N. 
Let us write down few first expressions: 




(3.13) 



(3.14) 




(3.15) 



19 



(3.16) 



Si = i(AWA[^U|l + ^fU|U[^l + 414141 - 41414=^1), (3.17) 

Trivially, the extreme values for Si and S2 are ±1. Then it is easy to argue from the construction of 

the sequence that also all mean values of S'^2) above 1 or below —1 cannot be explained with Local 
Realism. In deterministic theories based on LHV one of terms in l|3.15p vanishes, the magnitude 
of the other is cancelled to 1 with the fractional factor ^ appearing in each iteration. 

The quantum-mechanical bound was shown in [S] to be 2^^"^^/^. Thus so-called Mermin- 
Ardehali-Belinskii-Klyshko (MABK) inequalities turn out to be violated in some cases with a 
higher ratio than ever possible for Mermin or Ardehali inequalities. 

3.5 WWWZB Inequalities 



The derivations described above, especially due to Belinskii and Klyshko [5], utilize the com- 
plementarity mentioned above, however, only a certain configuration of such pairs of such comple- 
mentarjl^l expressions is used and each consideration leads to one certain inequality. 

The complete (in the sense explained further) set of Bell inequalities for N qubits given that 
each observer chooses between two dichotomic observables was derived by Werner and Wolf [88] . 
and independently by Weinfurter and Zukowski [87j and by Zukowski and Brukner |103j . They 
are known as WWWZB inequalities. 

The construction given in |103j is similar to the original argument in [6]. Let us consider 



all possible products (n^i(^i''' -^5^4'"')) with 



: ±1. Local- realistically, only one such 
For a particular choice of signs, s^*^! = 



expression at the time has non-zero value at the time 

/ 

is equal to ±2^. Thus the general Bell inequality for N qubits and 2 observables per site reads 



(A^i^) lhv{a12^)lhv (we remind that in these theories these mean values are also ±1), the product 



E 



N 



[fc] 



5W4'=i) 



< 2 



N 



(3.18) 



and is equivalent to 2^ inequalities 



sin 



.,s["l=±l 



N 



5W4'=l)) <2 



(3.19) 



where the sign function S{.) depends on signs sl'^l and takes values ±1. Within the set, we also 



find trivial inequalities, like 



ntr4 



[k] 



highest violation ratio is observed when |(3 



< 1. But also true non-trivial Bell expressions. The 
TqI reconstructs the MABK inequalities. 



It is possible to give a sufficient condition for an arbitrary state of N qubits to satisfy i|3.19p . 
Let us first rederive the necessary and sufficient condition for two qubits originally presented by 
the Horodecki Family [52], but rather in the formalism of [103) . The general two-qubit inequality 
can be written as 



2 

E 



((4i + (-i)''^'-4i)(4i + (-i) 



4i: 



< 4. 



(3.20) 



Recall that in Quantum Mechanics the mean values of our interest are computed from the corre- 
lation tensor, Ej^^„^j^ = (A'-^J...^':^1) ^ f ■ {d^l-^ <Si ...<Si dfJXji, ...Jn = 1,2)|!|. Both observers can 



^In Local Realism. 

^The scalar product of two tensors is here understood as A ■ _B = J^^^ 



20 



choose such local coordinate systems that the unit vectors defining their observables will satisfy 

^ [fe] , ^ [fe] „ ^ [fe] 

parallel to axes x, 
into ijOOl) we get 



a + d^2^ — 2e f^' cosafe/2 and a^^^ — a 2*^' = 26 2*^' sinafe/2, given that 61,6*2,63 are unit vectors 
parallel to axes x,y and z, respectively, and ak is the angle between a^^^ and al^^. Putting this 



(|Tll|,|Ti2|,|T2i|,|T22|)- 

(I cosai cosq;2|, I cosai sina2|, | sinai cosq;2|, | sinai sina2|) < 1, (3-21) 

where TyS are elements of the two-qubit correlation tensor (see Section 2.4). This means the 
following. One can take a correlation tensor of a state, initially given in any local bases. Then 
one constructs a modified tensor, T™°'^ in which all elements have been replaced by their moduli. 
If there are local coordinate systems in which any element of T™°'* exceeds 1, that is if T™°'' is 
not a correlation tensor of any physical state, the state violates p.20p . This greatly resembles 
Peres-Horodecki criterion [67l |49] (to be explained in Section 5.2), in which physicality of a state 
is also questioned after a certain operation, i. e., a partial transposition. 

The two vectors appearing in l|3.2ip can be made parallel by a proper choice of ai i a2 and 
coordinate systems. Moreover, the second vector built of trigonometric functions is always nor- 
malized. Hence, on the basis of the Cauchy inequality, \a-b\ < \a\\b\, the necessary and sufficient 
condition for a state to satisfy l|3.2ip reads 

2 

E ^afc < 1- (3.22) 

a, 6=1 

In the same fashion one can consider p.lSp for more qubits, up to the point in which we 
paralleHze the two vectors. This is in general not possible for more qubits due to the Hnear growth 
of free parameters and the exponential growth of of vector components. The analogous condition. 



E 



TLn < 1, (3.23) 



is only sufficient. If l|3.23p holds in all coordinate systems, the state can never violate any of 
inequalities l|3.19p . However, the condition with the modified condition, that is that 



max|r™°?|<l (3.24) 

holds in all local bases is still necessary and sufficient for p.lSp to be satisfied. The maximum is 
taken over all three-dimensional rotations, both before and after the modification of the tensor, for 
all local coordinate systems (each observer fixes his/her own Cartesian system). 

3.6 Drawbacks of WWWZB Inequalities [T04l [74] 

An example of a state, which neither satisfies l|3.23p . nor violates p.lSp is a so-called noisy W 
state, given by 

Pw,n{Vn) = Vn\Wn){Wn\ + (1 - ^iv) ^y , (3.25) 
|TyA,) = l(|10...0) + |01...0> + ... + 10...1)), 

which is certainly entangled for Vn > (^_i)^]v-i^.jv ' paper [72] shows that for this interval 
of Vn noisy W states reveal a violation of Local ReaHsm in a certain protocol. Namely, we 
demand that N — 2 observers perform 0-3 measurements on their qubits and obtain "+1" as their 
results. Such a result came with a high probability from the W admixture and with a relatively 
low likelihood from the white noise. Thus the overall W-state-to-noise ratio increases with every 
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projection. If this ratio is higher than after N — 2 projections, the state of remaining 2 qubits, 
on which no measurements were made (possessed in many copies) , violates a Bell inequality. Using 
this protocol we have shown that noisy W states lead to a stronger violation of Local Realism 
than noisy GHZ states. For iV > 11 we need a smaller admixture of the W state to white noise in 
order to violate a Bell inequality than a required amount of the GHZ state. This advantage grows 
exponentially with N. For example, for 3 qubits, the threshold for Vn is « 1/1.5523 w 

0.644212. 

More importantly, in [74] we show that the maximum of ^ ^^-^ T^bc ~ 1^3^ ^ Pw,z^z) or 
more generally, (3 — for an arbitrary value of TV. However, we have performed a numerical 

optimization to compare the value of V-i, above which (|3.18p is not satisfied to the one, above 
which (|2.15p is violated after the application of the protocol. We have shown that a pure 3- 

qubit W state violates l|3.18p by a factor of 1.523 < « 1.527. Thus there exists an interval 

0.6547 w < P3 < ~ 0-657, in which the condition ((3?23l) is not satisfied, but ifSTTSl is. 

Gisin [39] and Gisin and Peres [40] have shown that all pure entangled states violate the CHSH 
inequality. The general two-qubit state reads 

= coolOO) + coilOl) + ciollO) + ciilU). (3.26) 

The state is normalized to X^l j=o I'^ii P — 1- Now, we can find a normaHzed product state |0'0'), 
so that the scalar product (0'0'|?/') is real and maximal. This impHes that (l'0'|i/') = (0'l'|i/') = 0, 
where |1') are orthogonal to |0') in both local Hilbert spaces. Finally, since the global phases 
determining |1') are yet unset, we can choose (I'l'lV") to be real and positive. Thus we can rewrite 

as 

1^) =cosa|0'0') +sina|l'l'), (3.27) 

with < a < 7r/4. Such a procedure is called the Schmidt decomposition. It aims to minimize the 
number of non-zero coefficients of a state, or, more generally, a vector living in a tensor product 
of Hilbert spaces. It can be generalized to any bipartite system, or any number of qubits. 

It is simple to argue that the Bloch vectors of l|3.27p . (1?'^') and ((?'^^), have only 2:-components 
non- vanishing, T03 = T30 = cos^ a — sin^ a = cos 2a, and that the correlation tensor has only three 
non-zero elements: Tn = — T22 = 2 cos a sin a = sin2a,r33 — 1 (see Section 2.4 for the definition 
of Tijs). If we apply l|3.22p . we get Tf^ -|- T^^ = 1 + sin^ 2q! > 1 and = 1 only for the case of 
a product state, for which a = 0. By applying transformation reverse to the ones used in the 
Schmidt decomposition this proof can be generalized for all pure entangled states. 

Despite the result of Gisin and Peres, one can ask whether the WWWZB inequalities, despite 
of forming a complete set of inequalities for iV-partite correlation functions with two alternative 
observables per site, are violated by all pure states. Zukowski, Brukner, Laskowski, and Wiesniak 
citewl have shown, this is not the case. We have studied a family of states, which, despite of their 
obvious non-classicality, do not violate any of l|3.19p inequalities. The family are generalized GHZ 
states, 

\GHZ{a,N)) =cosa|0)^^+sina|l)^^ (3.28) 

with a as before. Scarani and Gisin [72] have shown that such states never violate MABK inequal- 
ities for sin 2a < 2^^"^'/^ with N > 2. This bound is vaHd also for WWWZB inequalities for odd 
N. The non-vanishing elements of the A^-particle correlation tensor are: 

T1...1 ~ —T221...1 = —T212...1 = ■■■ = 722221....! = ••• = sin2Q;, 

1 + f-l)^ 1-f-l)^ 
T3...3 = \ ^ + cos 2a. (3.29) 

The first Hue contains all the elements with an even number of subscripts 2, rather than 1. Thus 
the optimal sums in l|3.23p are max{l, 2^~^ sin^ 2a} for N odd and max{l -I- sin^ 2a, 2^~^ sin^ 2a} 
for N even. In the N odd case a possible choice is to consider the xz part of the correlation 
tensor, in which we have only two non-vanishing elements, r3...3 = cos2a,Ti...i — sin 2a, thus 
n=i 3^a...n — ^- Alternatively, we can consider the xy part of the tensor, which has 2^~^ 
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non-zero elements, all of modulo sin2Q!. Taking into account that 2^~^ sin^ 2a is less than 1 for 
sufficiently small a, we have shown that there exist pure entangled states of odd N qubits, which 
never violate any WWWZB inequalities. 

It is also worthy stressing that for N even the optimal sum is always above 1 (except for a 
trivial case of a product state). For the purposes of the demonstration of this statement, let us 
focus on the case of = 4 described in |1Q4) and the xz part of T, which reads as follows: 

T[xz] = e 3 (g) 6*3 (g) 63 (g) 63 + sin 2Q!eri (g e i (g e i (g e i . (3.30) 

First, let three of the observers rotate their coordinate systems by tt/A around the y-axis. Then 
the correlation subtensor takes a form of 

%.] - 2-3/2 ((g., ^ g., ) ^ (g., _^ g., ) ^ (g., ^ ) ^ 

+ sin2a(e i - e^,) g) {e[ - e^,) g) {e[ - e!,) g) e[) . (3.31) 
Next, we replace the elements of the tensor with their moduli: 

f^'i = 2-3/2 (^g^ + 6-3')® (e/ + 63') ® (e/ + 63') (g (sin 2aei' + 63') . (3.32) 

Finally, first three observers shall perform such transformations that — -^(e/ + 63') and the 
fourth observer should define e{' as , ^ (sin 2ae{ + eV). The tensor has now only one non- 

-^Z 1+sin^ 2a 

zero entry: 

T|'™f'* = Vl + sin2aei" (g e/' (g e/' (g e/'. (3.33) 

Applying l|3.24p to l|3.33p we see that at least one of WWWZB inequalities is violated by all such 
states for < a < 7r/4. The same argument can be applied for more qubits. The inequality, which 
is always violated a generalized CHSH inequality: 

IN N \ 

( (^f' + 4'') n ^i"' + - 4'') n ) ^ 2. (3.34) 

\ fe=2 fc=2 / 

The two diflFerent behaviors, for TV odd and even, are related to an interpretation of 
n=i '^a...n- Following [20j we argue that this sum can be understood as the amount of informa- 
tion stored in some of correlations between the qubits. This interpretation is justified by the fact 
that when we measure an observable A with two outcomes, "±1", which have equal degeneracies, 
1 — A2(A) — (A) 2 is equal to 1 in the case of one of the results occurring deterministically, and 
when both results occur with equal probabilities. The first situation would suggest that we have 
learned one bit of information about the systerr0- Product states, and hence their statistical mix- 
tures, cannot contain more than one bit encoded in these correlations. Thus if the condition p.23p 
is not met, the presence of entanglement is immediately impHed. In such a case, it is possible, but, 
as we have mentioned, not necessary, to violate Bell inequalities. 

3.7 WZLPZB Inequalities 



A natural way of extending the set of known Bell expressions is to consider more than two 
observables per site. An example of such an inequality was given for three qubits by Wu and 
Zong in [96]. The inequality given therein utilizes four alternative measurements for two of three 
observers and two measurements of the third one. A generalization of the inequality from |96] 
was demonstrated by Laskowski, Paterek, Zukowski, and Brukner [59j. They have noticed that for 
LHV models both the expressions, 

^12,5,- E ^i(fci,^i)(4'' + (-i)'^4'')(4" + (-i)'^4") (3.35) 

/ci,ii=l,2 

'^We need to stress, however, that it is crucial that the degeneracies are equal. Consider an extreme case, in which 
"+1" appears much more often in the spectrum of A than " — 1". For a maximally mixed state, which obviously 
contains no information, (A) would then be close to 1. 



23 



and 



^34,5.- E ^2(fc2,;2)(4"+(-l)'^4')(4"+(-l)'^4"), 



(3.36) 



fc2,i2 = l:2 



are always equal to ±4, similarly to ^ being equal to ±1. Thus a new three-qubit inequality can 
be derived in a following way: 



|(^12;12,34)| 



< 16. 



(3.37) 



Since any of the sign functions Si, 82,8 can have one of 2* = 16 forms, l(337)) stands for (2^) 2^2 
inequalities. If none of S, 81,82 is factorisable with respect to (— I)'', (— 1)',; (— 1)*^% (— 1)'% or 
(— 1)'^^, (—1)'^, respectively, the inequalities are equivalent to 



(AW(^fl + 4^1) + 41(Afl - 4'l))(Afl + 4^1)+ 



(41(4^1 + 4^1) + 41(4^' - 4^'))(4^' - 4^') 



< 16, 



(3.38) 



which was found in [96j. If or 82 is the only factorisable function under the consideration, we 

get an inequality equivalent to p.38p with A^^^ — ^2^' or A\^^ — A\1'^ . In other cases, as well as if, 

for example, A^^' = Ag^' , ^2^' = 4^' ) ' = ^3^' > ^2^ = ^4 ' ■• the expression already belongs to the 
set ((3A9l) . 

It is now possible to extend the derivation to more qubits, for examlpe, for = 4: 



l(^12;;12; 



12,34;:34;56J8/ 



EL=l^(fc,0(^12;12,34 + (-l)'=A34;56,78)(4"+(-l)'4") 



< 16. 



(3.39) 



We thus obtain a series of A^-qubit Bell inequalities, in which the last observer chooses from 2 
apparatus settings, and the previous has twice as many as the next one, except the first two, who 
both have a choice of 2^~^ settings. 

We recall that for WWWZB inequalities we chose the sum and the difference of the two vectors 
defining a pair of observables for each observer to be proportional to ei an 62 I'^' . The same 
can be now done for the vectors of the last observer. The mean value of the Bell operator is then 
expressed as two terms, dependent on disjoint set of observables of all other observers. Thus such a 
procedure can be applied for the next observer, with independently chosen local Cartesian systems 
in each term. Repeating the reasoning from Section 3.5 we finally reach the necessary and sufficient 
condition for a state to satisfy the inequality. For the three-qubit inequality l|3.38p it reads 



2 

Tabl 

a,b=l 



< 1. 



(3.40) 



,6=1 



Carol's coordinate system is the same in both therms. Cartesian bases of Alice and Bob for the 
second term may be chosen different from those for the first one. The analogous condition for a 
state to satisfy l|3.39p reads 



a,&=l,2 



rp2 I \ " rpl2 I \ " rpll2 , \ " rplH2 ^ 1 



(3.41) 



a,fc=l,2 



a, 6=1, 2 



a, 6=1, 2 



In any of the four terms the two first experimenters are allowed to arbitrarily choose their coordinate 
systems, individually for each term. The third observer has one Cartesian basis for the first two 
terms and some other for the other two. The forth works in a fixed coordinate system in all four 
terms. 

Importantly, conditions l|3.40p . p.4ip . and similar ones we would obtain for more qubits are 
necessary and sufficient. Thus the WZLPZB inequalities are violated by noisy W states p.25p for 

any Vn > l/i/3 — unlike in the case of WWWZB inequalities. 
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All generalized GHZ states considered in |104j violate the WZLPZB inequalities for a 7^ 0, 
also for TV odd. To show that, let the observables of the last observer be cos^iCTg^' ± sin6'i(T[^'. 
Now, only such elements of the iV-particle correlation tensor enter the condition for satisfying or 
violating the inequality, which have the last subscript equal to 1 or 3. In the first case all others 
should choose observables in such a way that the condition is entered by the elements of the tensor, 
for which the subscript "2" appears an even number of times, while the remaining subscripts are 
"1", e.g. T11...1 or T2211...1. Since the value of the last index is fixed, there is 2^~^ such terms, each 
equal to ± sin 2a. The total contribution of these elements to the sum is at most 2^-2 sin^ 2a. If 
the last subscript is equal to 3, the rest of observers should ensure that 3 = cos^ 2a enters the 
condition. The total sum would then be 1 + (2^~^ — 1) sin^ 2a > 1. 



3.8 Convex Hull Problem 



As first considered by Froissart [36], Bell inequalities define hyperplanes in a statistical spac^. 
In between all these hyperplanes for a given experimental setup there is a convex polytopqj. This 
polytope contains all possible statistical distributions which can be explained with LHVs. This 
makes reasonable to ask two following questions about Bell inequalities. 

The first one is "Is the family, to which a Bell inequality belongs, complete?". That is, if the 
family completely bounds the polytope of LHV statistics in its most compact form. The answer 
to this question depends on the dimensionality of the statistical space. For example, WWWZB 
inequalities are complete when each observer can choose between only two measurements, but, as 
we have shown, there are physical probability distributions, which satisfy the inequalities (in this 
case, belong to the polytope), but violate WZLPZB inequalities. 

The other problem that can be addressed is whether or not a given inequality is tightly. 
The original Bell theorem (see Section 2.4) states only that there exist experiments, quantum- 
mechanical results of which cannot be explained within Local Realism. For the purposes of the 
falsification of LHV-based theories it is enough to find any hyperplane, which separate any phys- 
ically accessible point of the statistical space from the polytope of statistics explainable by LHV. 
However, for the sake of detecting of useful entanglement in as many states as possible, we need 
to find optimal, or in other words, tight Bell inequalities. The conditions for tightness are that the 
hyperplane contains at least D extreme points of the polytope and that the whole its interior lies 
at the same side of the hyperplane. D stands for the dimensionality of the space. This problem is 
called convex hull problem. It is thus the usual procedure to find tight Bell inequalities: to find a 
set of the extreme points of the LHV-permitted polytope, and then to list hyperplanes containing 
of D of them, which do not intersect the interior. Such lists were given, for instance, by Pitowsky 
and Svozil [59], Sliwa [78], and Collins et al. [13 . It also turns out that [S [13 [HI [3 [lOHl EH] are 
also elegant, yet unconscious, derivations of inequalities optimal in this sense. 

Out of many examples from [69], let us give the one for two qubits and two events per site, 
denoted as, Ai, A2, Bi, B2. A point of the statistical space is described by 8 probabilities arranged 
to a vector: 

(P(Ai), P(A2), P(Bi), P(B2), P(Ai, Pi), P(Ai, P2), P(A2, Pi), P(A2, P2)). 

*By a statistical space we shall mean a real vector space with components given by probabilities of certain events 
(then the space is referred to as a probability space) or mean values of certain operators (as we will consider products 
of local observables, this will be called a correlation space or a mean values space). 

polytope is a finite region of a hyperspace bounded by a finite number of faces. The word "polytope" closes 
a sequence of notions "point, line segment, polygon, polyhedron,...". 

^"A tight inequality defines a hyperplane, which actually contains a face of the polytope. 
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Local realistic theories describe sixteen points of this space: 

(1,1,1,1,1,1,1,1), (1,1,1,0,1,0,1,0), 

(1,1,0,1,0,1,0,1), (1,1,0,0,0,0,0,0), 

(1,0,1,1,1,1,0,0), (1,0,1,0,1,0,0,0), 

(1,0,0,1,0,1,0,0), (1,0,0,0,0,0,0,0), 

(0,1,1,1,0,0,1,1), (0,1,1,0,0,0,1,0), 

(0, 1, 0, 1, 0, 0, 0, 1), (0, 1, 0, 0, 0, 0, 0, 0), 

(0, 0, 1, 1, 0, 0, 0, 0), (0, 0, 1, 0, 0, 0, 0, 0), 

(0, 0, 0, 1, 0, 0, 0, 0), (0, 0, 0, 0, 0, 0, 0, 0), 

and the complete set of inequalities reads 





< P{A, 


B,)<P{A,),P{B,)<1, 




(3.42) 


1 < P{AiBi) - 


f P(AiB2) + 


P{A2Bi) 


- P{A2B2) 


-P{Ai) 


- PiBi) < 0, 


(3.43) 


1 < P{AiB2) - 


f P(AiBi) + 


P{A2B2) 


-P{A2Bi) 


-PiAi) 


- P{B2) < 0, 


(3.44) 


1 < P(^2Si) - 


f P{A2B2) + 


P{AiBi) 


-P{AiB2) 


- P{A2) 


- P{Bi) < 0, 


(3.45) 


1 < P{A2B2) - 


f P(A2Bi) + 


P{AiB2) 


-P{AiBi) 


- P{A2) 


- PiB2) < 0, 


(3.46) 



i,j — 1,2. The necessary and sufficient criterion for a point to be in agreement with Local Real- 
ism is that it satisfies all these inequalities. (|3.42p are trivial inequalities, which define properties 
of probabilities, whereas l|3.43p . I|3.44p . p.45p . and p.46p are Clauser-Horne inequalities, equiva- 
lent to (|2.25p . However, permitting observers to consider three events leads to more compHcated 
inequalities given in [69]. 

3.9 Tight Bell Inequalities with up to Three Settings per Site 



It is important to stress that the approach of convex hull to the problem of finding Bell inequal- 
ities is suitable for both the probability space and the correlation space. Derivations of Pitowsky 
and Svozil [69] were made for probabilities, whereas Zukowski [T02j attempted to find tight in- 
equalities by a direct analysis of the polytope spanned in the correlation space. The subject of his 
interest was a situation, in which each of iV observers chooses between three dichotomic observables 
with usual outcomes, " ± 1" . 

For the sake of illustrating the problem Zukowski considers the case of iV = 1 thus working in 
the space of vectors built of mean values of observables vIq^', ^f', and ^2^'. The inequalities found 
for iV = 1 cannot be interpreted as Bell inequalities, since tho whole physically allowed fragment 
of the space is also per definition accessible with local hidden variables. 

An elementary event is that measurements of the three observables give results Oq', a^^', Oj^^, 
respectively. The probability of such an event shall be denoted as P{a\l\ a^i \ a^2^). The mean value 
vector shall be defined as 

E= y: (4^Ul^i,4^V(4^',4^',4^'). (3.47) 

The polytope is a cube, faces of which connect the following quadruples of points {a,b — ±1): 

(l,a,6), (-l,a,6), (a, 1,5), 
(a, -1,6), (a, 6,1), (a, 6,-1). 
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For a given face any three vertices form a complete linearly independent set in 3?^, however, vectors 
of this set are not orthogonal. Any other vector belonging to the space can be uniquely expressed 
as a linear combination of the three. In particular, the fourth vector of the face can be decomposed 
to the first three with coefficients of moduli 1 and which add up to 1. For example, 

ab{l,a,b) = (0,0,0). (3.48) 

a,b=±l 

More generally, we recall the fact from Analytic Geometry that a hyperplane in a ZJ-dimensional 
space which does not cross the origin of the coordinate system, is defined by D linearly independent 
vectors. All other vectors which point the hyperplane are such linear combinations of the initial 
vectors, that the expansion coefficients add up to 1. Rephrasing this fact explicitly, we take a 
complete set of D lineary independent vectors, {rijlLi. ^ hyperplane is a defined as a set of all 

such vectors r = J2iLi ^i^i^ that J2iLi ti = ^■ 

Now, to get more intuition on the problem for more complicated cases, Zukowski suggests to 
consider three randomly chosen vertices of the cube. If these three vertices belong to the same face 
of the cube, the fourth vertex at the plane can be expressed as a linear combination of the three 
with coefficients adding up to 1. 

In the other case, the three vectors are linearly independent, but do not form a face. Thus 
there must exist a pair, v and —v, such that its elements lie at two sides of the hyperplane. The 
hyperplane cuts the polytope through. 

The last possibility is that the three vectors do not estabfish a complete set. Thus there are 
vector V and —v which are not expressible by the trio. Moreover, such a trio must consist of a pair 
v' and —v', thus the hyperplane also cuts the cube through. 

Finally, Zukowski presents the following corollary. If one builds a non-orthonormalized basis 
out of three vectors representing the vertices of the polytope, which defines its face, all other 
vertices are described by convex combinations of basis vectors with coefficients equal to ±1. 

Let us define the analogous problem for two qubits. The correlation space is described by a 
correlation matrix, 

i: E (4^' , , 4^' ) ^ (a^^' , , (4^' , 4^' , 4^' ; a^^' , , ) ■ 

(3.49) 

We thus work in a nine-dimensional correlation space and the polytope is spanned between 2^ — 32 

verticeini- We can always divide all the vectors representing them into two groups in such a way 
that one element out of each pair {v \ —v''} belongs to one group, and the other to the second. In 
each such group out of 16 vectors we can choose 9 vector forming a complete set. Which such sets 
then define a face of the polytope? 

Zukowski now notices Fact 1. He considers a randomly chosen complete basis of vectors {v^}^^^ 
and takes an assumption that there exists seven more vectors, which have the expansion coefficients 
summing up to 1. Thus these sixteen vectors constitute one hyperplane. Zukowski now aims to 
prove that no other vertex belongs to the hyperplane and all of them lie at the same side of it. 

Proof: The necessary and sufficient condition for a vector x to be part of the corresponding 
hyperplane is that the quantity 



D[x,v^, ■.■,v^] = det 



is equal to 0. Moreover, D[x, v^, v^] is always positive at one side of the hyperplane an negative 
the other. 

^^Having in total six measurements, one might expect to have 2^ = 64 vertices. Please note, however, that 

{a]^^ , aj"""' , aj^' ) ® (<i[f' , , Hj^' ) = {— Oq^' , —o-^i^ , — 02^' ) ^ (^i^cf' 1 ^o^i^' 1 ~'^^2^)j cHid thus their number is decreased by 
a factor 2. 



Xi X2 



vl 



(3.50) 
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The vertex vq, which is from outside the hyperplane is either given by —v''{k = 1,...,9) or 
tkv'^ with X]fe=i — —1- In either case, we can add such a linear combinations of rows 2 
to 10 to the first row that only its first entry is 2, while the rest of them vanishes. It is then clear 
that 

D[v°,v\...,v'^] = 2d[v\v^,...,v'^] (3.51) 



with 

d[v^, ...jV^] det 



vl ... 



(3.52) 



QED. 

Fact 2 states that plane containing both v and —v contains also — (v + (— 1/))/2. 

Fact 3 is that coefficients of expansion of any vector in a non-orthonormal basis are integers. 
In fact, we can show that in our case these coefficients are ±1 or 0. 

Proof: let us choose 9 basis vectors as (l,a',6') (8) (l,c',(i') where a',b',c',d' = ±1, but never 
{a',b') = (1, 1) or {c',d') = (1, 1). Now, we use l|3.48p to obtain 6 more vertices: 

(1,1,1)® (l,c',d') = {l,a',b')(g)(l,c\d'), (3.53) 

(a',b')#(l,l) 

(1,1, 1)^(1, 6', a')- (l,a',6')® (l,c',d'), (3.54) 

and the 16th one: 

(1,1,1)0(1,1,1)= Y (l,a',6')0(l,c',d'). (3.55) 

{a',b'),{c',d')^(l,l) 

Sixteen other vectors are obtained by a global sign flip. For a set of nine in which other pairs of 
(a', b') or (c', d') being excluded, a similar argument can be presented. If a choice of nine is more 
complicated, yet the set is complete, it is sufficient to use transformations inverse to p.53p . I|3.54p . 
and (|3.55p . coefficients of which are also ±1 or 0. Applying these transformations once more, we 
see that all coefficients must be integers. QED. 

Finally, assume that we have found a hyperplane containing between 9 and 15 veritces, i. e., 
there exists a pair of vectors, v " and — u neither of which ends at the hyperplane. Fact 4 is that 
they lie at two opposite sides of the hyperplane and the latter cannot define a face of a polytope. 

Proof: obviously, we can expand ztiP — ± ^kV^- By adding an appropriate combination 

to the first row of the matrix used to compute l|3.50p we get 

9 

D[±v'',v\...,v^] = lTY^kd[v\...,v% (3.56) 

fc=i 

The previous proof showed that all coefficients are integers. It is, moreover, easy to argue that 

they must add up to an odd integer. Since it is impossible that J2t=i^k — ^, D has different 

signs for v° and — u °. The hyperplane must cut the interior through. We hence conclude that all 
faces of the polytope contain exactly 16 vertices. QED. 
The Bell inequalities will be given by 

±D[E,v\...,v^] <0, (3.57) 

where the global sign is to be determinecf^. and are 9 of 16 vectors defining a hyperplane, 

which has the properties emerging from the described facts and contains a face of the polytope. 

Now, we are only left with the problem of finding valid sets of vertices. One such set is certainly 
{(1, a, b)(^{l, c, d)}a,b,c,d=±i- Other sets, claims Zukowski, are given by S{a, b, c, d)(l, a, 6)(Xi(l, c, d), 



^^Certainly, the maximally mixed state, ^l23x23i which carries no information, satisfies all Bell inequalities. 
The global sing of the left-hand side of H3.57|l shall be hence chosen the opposite as of d[v^ , ...,v^]. 
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with a sign function S{a,b,c,d) such that X^a h ^^"^(^^^fi) = ^cdS{a,b,c,d) = 0. This also 
guaranties that, e.g., J2abc o,bcS{a, b, c, d) = J2a bed <ibcdS{a, b, c, d) = 0. The exphcit form of the 
sign function is 

S{a, b,c,d) ^ X + Aa + Bb + Cc + Dd + Eac + Fad + Gbc + Hbd, (3.58) 

where X, ...,H are constants. 

Proof: let us choose the basis as S{a' ,b' , c' , d){l, a' ,b') ® (l,c',d'), with cases {a',b') = (1,1) 
and {c',d') — (1,1) excluded. 

For the beginning, let us consider a constant sign function. The remaining 7 vertices belonging 
to the face are obtained from l|3.48p . and therefore the coefficients of expanding (1, 1, 1) (g) (1, c', d') 
and (l,a',6') (g) (1,1,1) add up to -J2{a'b')^(ii)a'b' = - E(c' i) = ^- Eventually, 
(1, 1, 1) ® (1, 1, 1) = E(c',d')#(i.i) a'b'c'd'{l,a', b') ® (1, c', d') = 1. 

We will now generalize these considerations for an arbitrary sign function. When S{a, b, c, d) is not 
constant, we cancel its action by putting it into the coefficients: 

S{1, 1, c', d'){S{l, 1, c', d')(l, 1, 1) ® (1, c', d')) 

= E(a',6')#(i,i) a'b'Sia', b', c', d'){S{a\ 6', c'd'){l, a',b') ® (1, c', d'), (3.59) 
5(a',5',l,l)(5(a',6',l,l)(l,a',6')® (1,1,1)) 

= E(c',<i')#(i4) c'd'S{a\ b', c', d'){S{a', b' , c'd'){l, a', b') ® (1, c', d'), (3.60) 
5(1, 1,1,1)(5(1,1,1,1)(1,1, 1)0(1,1,1)) 

= E(a',ft')#(i,i) E(c',d')^(i,i) <^'b'c'd'S{a', b', c', d')(^(fl', 6', c'd')(l, fe') ® (1, c', d')- 

(3.61) 

Thus the sums of the coefficients of the expansion are given by 

^S{\, 1, c', d') E(a',6')#(i,i) <^'b'S{a', b', c', d') = S{1, 1, c', d')' = 1, (3.62) 
-5(a', 6', 1, 1) E(c',d')#(i4) c'rf'5(a', b' , c', d') = 5(a', b', c, df = 1, (3.63) 
-5(1, 1, 1, 1) E(a',6')#(i,i) E(c',<i')#(i4) «'^Vd'5(a', 6', c', d') = 5(1, 1, 1, If = 1, 

(3.64) 

where these identities follow from l|3.48p and (|3.58p . 

Zukowski next considers a sign function cr(a', 6', c', d') which does not have the form of i|3.58p . 
In particular, let us focus on the case of E(c' d')^{i i) c' d' a{a' ,b' ,c' ,d') = xa{a' ,b' ,1,1), x ^ —1. 
This means that either a; = 1 or x = ±3. As the relation 

E(c'd')#(i4) ^''^''^(°'' '^^ d'){o{a', b', c', d'){l, a', b') ® (1, c', d')) 

= -a{a', b', 1, l){a{a', b' , 1, 1)(1, a', 5') ® (1, 1, 1)) (3.65) 

holds, the sum of the expansion coefficients is 

^ c'd'a{a',b',d,d')cT{l,l,c',d') = -xa{a',b',l,lf ^^x^l, (3.66) 

i. e., the hyperplane does not contain the vertex. QED. 

The Bell inequalities in the form of l|3.57p look complicated, not suitable for direct applications. 
The equivalent form for the case of two qubits and three observables per qubit is 

S{a,b,c,d)E-{l,a,b)(g>{l,c,d)<2^, (3.67) 

a,b,c,d—:tl 
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or, more explicitly, 

Ea,6,c,d=±i'5'(a,^c, d) 
X {Eqo + aEiQ + bE2o + cEqi + dEo2 + acEn + adE^ + bcE2i + bdE22) 

< 24, (3.68) 

or, using l|3.58p . 

-'^^'00 + CEqi + DEq2 

+ AEit^ + EEii + FEi2 

+ BE2Q + GE21 + HE22 < 1, (3.69) 

where = (^fUf ). 

Such a form of the inequality can be straight-forwardly generalized for more complicated cases, 
hke N qubits, or more observables per site. 

3.10 Explicit Form of Inequalities |93] 

Zukowski |102) has presented a general method to find tight Bell inequalities, but never gave 
the inequalities in their explicit form. The problem was thus brought under further analysis by 
Wiesniak, Badzi^g, and Zukowski in [93j. The paper aims to show a method of building a sign 
function. 

We shall start with the case of A'' = 2 For completeness of the analysis we first need to redefine 
the sign function. In the discussed case, instead of two, let S be dependent on three signs per 
observer, so that p.68p has the form 

^ E E Si4\s^\4^;s^^\s^\s^^^)(±s^l^sfE^]<l, (3.70) 



U ..,[11 ..,111 =± 1 ,421 ,4^1 =± 1 \ i J=0 



and the discrete Fourier transform (up to a multiplicative constant) of the sign function is equal 
to 

Sis^\^,^;s^,^,^)^±g.,s^^. (3.71) 
Coefficients Qij are equal to 

9.-^ E E Si^,^,^;^,^\^)^sf (3.72) 



11 Jl] Jll— n J2] J2] 12] 



Sq '^1 '^2 



so that p.70p can by shortly written as 

2 

E 9^JE^J < 1. (3.73) 

Let us concentrate on the properties of the coefficients. We start with defining functions, which 
we shall call deltas of some order. First order delta with respect to s[^' is given by 

A ™ = (3.74) 

Since for any given combination of ah other signs the sign function may, or may not flip, A^[i] 
takes values of ±1 or 0. Similarly, we define second order delta with respect to and s^^': 

A^ui(4'l = l)-A^„(4" = -1) 
Aju j.i = ^^^^ . (3.75) 
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Under the flip of Sj A^m can remain unchanged, or change between ±1 and or between 1 and 

— 1. Thus possible values of a second-order deltas are 0,±i,±l. In general fcth order deltas can 
take values between —1 and 1, which are multiples of 2^"'"'. 

The other obvious thing about deltas is that for TV qubits TVth order deltas are coefRcients of 
the Fourier transform of the sign function. The higher order deltas can be thus reconstructed by 



[2] 

■ 9ioSo 



n .[2] , „ J2] 



and 



[1] 



,[2] 



, O]^ 



A [1, 



A 111 4 



[1] 



Other facts, which would be helpful in constructing valid sign functions, are 



(3.76) 
(3.77) 

(3.78) 



The first equality follows from Fact 4 from |102j that each face of the polytope contains 16 ver- 
tices. The second is a consequence of Parseval's theorem, which states that the Fourier transform 
preserves the EucHdean norm of the transformed object. Since the norm of any sign function is 



E 



(3.79) 



also its norm in the Fourier picture must be constant. We have dropped the factor guarantying 
normalization of the Fourier transform, however. 

It is easy to confirm that the only non-trivial possibility for first order delta is ^(sf ± sf}) 
with i ^ i' . For simplicity let us take i = and i' = 1. Any more terms would allow the delta to 
exceed 1 for a certain choice of signs. Moreover, the norm of such a delta |A [i] p, that is the sum 

of squares of all coefficients entering the delta, is ^ . Thus from l|3.78p we conclude that two such 
deltas establish a sign function, one with a " -|- " sing, the other with " — " . The only non-trivial 
sign function for two qubits reads 



(3.80) 



which represents the CHSH inequality, (|2.15p . Thus we conclude that for correlations of two qubits 



the polytope is bounded by trivial inequalities 



< 1 and CHSH inequalities. 



Let us apply the same method for the case of three qubits and three observables per site. 
Coefficients gijk satisfy conditions similar to l|3.78p and take values 0, ±i, ±i, ±|, ±1. The last 
case, of course, corresponds to trivial inequalities and thus it will be skipped in further discussions. 
The Bell inequalities take a form 

ijk 



(3.81) 



where £;,,fe = (AfUf^Jf^). 

The following table gives all possible forms of A [i] [2]. It also gives a list of other A [i] [21 s, a 

certain one can go with in order to create a vaHd delta of the first order (obvious repetitions are 
avoided) : 
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• _2_ 



i(3(-irsi=^' + (-i)"4'') 



goes with: 



i(2(-l)",sP 

+(-irs|fi + (-i)°si'i) 



±i((-irs"F-(-ir4^') 



13K 



or 



i((-ir4=*j + (-ir4'0 



or 



and 

±i((-irsr-(-i)"4i) 



Tar 



alone or with 
or 

two ±i((-irsL'' - (-l)"s{f' )s 



i((-irs"F+(-i)"4') 



I3K 



any one of 
or 

±i((-irsL'i + (-i)»4J) 

and 



or 

J3] 



±i((-ir4' + (-i)°s?i) 

and 

±i((-irslfi-(-irs?i) 



±i((-irs"f +(-i)"4J) 



or 

,[3] 



and 

±3((-l)'"si" + (-l)°si") 



In the table (x, y, z) is a permutation of (0, 1, 2) and m, n, o = ±1. 

Thus any first order delta belongs (after local transformations, i. e., permutations of observables 



and sign flips) to one of families listed below (we take Sg ' — Sg^' 



J3] 



1 for a time): 
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A^[i] 


\A [111" 

■• 


An 




16 
Ifi 


At 


4 I O + *i + + )) 


12 
Ifi 


Att 


2 ' *1 ) 


8 
Ifi 


A m 




8 
Ifi 


A ri^ 




8 
Ifi 




1(2 + 5^3] ,^13] ,12]/ 131 _ 131 ^^ 


8 
Ifi 






4 
Ifi 


At/ I I 
^VII 


4VI>''1 ^2 J *1 \^ °l 7 ^2 ^-'^ ^2 I) 


6 
Ifi 


Ay III 


4 ''1 ^ ''2 ^ '^l V-^ ''1 / ''2 °2 1 1 


10 
Ifi 


Aix 




12 
Ifi 


Ax 


i(2 + 24^' + (.f'+4^')(l-.f')) 


12 
Ifi 


Axi 


i(2 + ,sf'+4^'+4^'(2-.f'-4'')) 


12 
Ifi 


Axil 


i(2 + 4^'(4^' + 4^') + 4^'(4^'-4'')) 


8 
Ifi 



Now, it is a necessary but not a sufficient condition that such a set of first order deltas enter 
the sign function, norms of of which add up to 1. Note that Aq is already a sign function, which 
represents a CHSH-Hke inequality: 

(-E'ooo + -E-ooi + £-010 ~ -Eon) < 2. (3.82) 

In other cases we need perform such local actions (observable permutations and sign fiips) on 
A [1] , A [1] , and possibly A [ij , that together they constitute a sign function. 

As we have mentioned, WWWZB inequalities, i. e. those, which utilize no more than two 
observables per site, are tight. Thus it should be possible to obtain them from the analysis of the 
sign function. For example, by using A/ and Ajy we can obtain the following sign function: 

+ s^^\l + sf^){l + sf^)), (3.83) 
and by putting it into p.8ip we have 

-(— 3i?ooo + £^001 + £^010 + £^011 

+Ewa + Ewi + Eiw + Eui) < 1. (3.84) 

Taking two A//s or two A/// we construct 

s = ln + sf^ 

+ 4\l~sf^)), (3.85) 

s = ln + sf^ 

+ 4"(l-4"), (3.86) 
which lead to CHSH-like inequalities l|3.82p . Another possibility given by these deltas is 

S - \il + sf 

+ 411421(1-43]), (3.87) 
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which implies 



2 (-^000 + -£'001 



+Eno-Ein) < 1. (3.88) 
A combination of A// and A/// leads to one of Mermin inequalities: 

s= + 

+ 4"(4"-4")), (3.89) 
^(-E'ooo + Equ 

+Eioi - Eno) < 1. (3.90) 

This closes the set of WWWZB inequalities, at least one of each kind, for three qubits. Now we 
ready to show Bell inequalities for which at least one observer has a choice of three observables. 
For example, two Ays can be arranged into 

s = l(2 + .ri + .r + .ri(4^'-4^') 

+ 4^1(2 -4^1 -4^1 + 4^1(4^] -4^1))), (3.91) 

what expresses the inequality 

^ (2£^ooo + -E'ooi + £-002 + £-011 — £-012 

+ 2Eif)o — Eioi — Eif)2 — Em + E112) < 1. (3.92) 

This sign function can also be a combination of two Axis or Ay/ and two A/xs. Also two Axiis 
can be used to construct 

s = i(2 + 4^1(4^1 + 4^1) + 4^i(4''-4'') 

+ 4^1(2 -4^1(4^1 + 4^1)- 4^1(4^' -4^'))- (3.93) 

which can also be obtained from A/77 and two Ay/s. After putting into (|3.8ip it gives 

^ (2i?000 + £'011 + £'012 + £-021 — £-022 

+ 2£^ioo — £111 — £112 — £121 + £122) < 1- (3.94) 
Let us also consider a pair of Ay// and A^: 

s ^ i(2 + 4n4^i + 4'i(i-4'') + 4^ki-4^') 

+ .}((4^i- 4^1) + 4^1(1 -4^1) -4^1(1 -4^'))). (3.95) 

With this sign function, which can also be obtained from Ay and two Ay/s, we can construct an 
inequality of a form 

^ (2£ooo + £001 + £002 + £010 — £011 + £020 — £-022 

+ £101 — £102 + £110 — £111 — £120 + £122) < 1. (3.96) 
The last possibility of constructing a sign function is with two Ay//s and a Ay: 

s = i(i + 4^i + 4^i(i + 4'') + 4^'(4^'-4^') 

+ 4^i(i-4'' + 4^'(i-4^')-4^k4''-4^')) 

+ 4^i(4^' + 4^'-4^'(4'' + 4^'))), (3.97) 
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which leads to 



^ {Eooo + £-001 + £-010 + E012 + E021 — E022 
+ EiQo — EiQ2 + Eiio — Em — £'121 + E122 

+ -E'201 + £'202 ~ £211 — E222) < 1- (3.98) 

Interestingly, (|3^ . ijaMj) . and ([3^ are special forms of (p^S)) derived in [96]. Namely by 
choosing ^[f' = A^^^ we obtain l|3.94p and subsequently l|3.92p with ^[f' — J^2^ . Finally, we obtain 
ijOe]) by choosing ^{f ' = A^^' and ^[f' = ' . Further such simplifications lead to WWWZB 
inequalities. 

ijaM)) . ijaMl) . and l(3:96| are also special cases of the 3 x 3 x 3 inequality, (|3^ . faMl) is 
obtained by e.g. A^^"^ = ±-4|,^'- With e.g. = iA^^ we obtain ij^Ml) . By putting e.g. 

A^^' = ±a[,^1 we get ((3:961) . 

It is now interesting to find conditions on states to satisfy these inequalities. This can be, of 
course, done in the same fashion as for WWWZB and WZLPZB inequalities and basing on the fact 
explained in the previous paragraph. The sum and the difference of two normalized vectors are a 
pair of observables are two orthogonal vector. It important to stress that a fixed direction of one of 
the original vectors does not fully determine the plane spanned by the pair. However, as observers 
choose fewer observables than in (|3.38p . the conditions presented below are only sufficient, but not 
necessary. Repeating the argumentation we obtain that l|3.40p is a sufficient condition for p.96p 
to be satisfied. 

Now, if one of the sign functions 81,82 in p.37p is factorisable, for example S'i(fci,^i) with 
respect to (—1)'^^ and (— we get i|3.94p . In such a case the condition reads 

Tin + E (^2y)' < 1, (3.99) 

where, again. Bob and CharHe can choose different coordinate systems in both terms. 

As for (|3T92l) . we can take ao^^' + ai^^' = 2 cos aei^^] , ajil - ail^l = 2 sin 06*2^ , ai'^l + a-^^'^ = 
2cos7ei['^l and a^^^ — a^'^ = 2sin7e2''^l (ei = (1,0,0), 62 — (0,1,0)). After performing necessary 
calculations similar to ones from Section 3.5 we obtain the first form of the condition: 

(7^S)' + r|n + (rii'^)'<l- (3.100) 

Superscripts (c) and (6) denote that in these terms Carlie and Bob, respectively, are allowed to 
perform arbitrary rotations of coordinate systems in which the second term is expressed. Please 
note that l|3.100p contains only 3 elements of the correlation tensor, and thus (|3.92p can be violated 
at most by a factor -\/3. The three terms can be brought to a common basis, since max(c) (T^^ii)^ = 
T^ni + Tlx2 + 7^n3> and similarly, max(b) (r2^i'^)2 = T^^^ + Tl^^ + Tl^^-. 

3 

r2n + E(^n^+^2'2) < 1. (3.101) 

i=l 

3.11 Additional Constraint on Local Realism |64| 

As we have mentioned in Chapter 2 the theories that are excluded by Quantum Mechanics [6] 
are based on two main assumptions; Realism, which states that physical systems posses properties, 
regardless if already measured, and Locality, which forbids superluminal effects. Nagata, Laskowski, 
Wiesniak, and Zukowski [64] have considered a case, in which a local realist is allowed to construct 
his theories only under a condition that the correlation function is a rotationally invariant function 
of unit vectors defining local measurements, i.e., 

E{a^,a^, a^)^f- {3^ ® ® ... ® a^). (3.102) 
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This form is independent of particular choices of local coordinate systems, but dependent on mutual 
relations between the measuring apparata and the physical system 

Let us first focus on equatorial observables, i. e., those, vectors of which lie in the xy-plane. 
The correlation function is now a function of angles, which determine positions of the vectors and 
can be seen as a vector living in the 3?^ space of real square-integrable functions. The measure is 
then given by dSl = 11^=1 ^"^i' the scalar product be 

{G,H)^ d4)i... ...d<j)NG{a^,...,a'^)H{a^,...,a'^). (3.103) 



We can now utilize facts from Analytical Geometry, especially the one that if for two vectors G, H 
there exist a vector F, such that 

(F,G) > (F,i?), (3.104) 

it is certain that G ^ H. Now, a Bell inequality can be constructed hj finding such a functional, for 
which, for a given state, {F,Eqm) > {F,Elhv), where Eqm and Elhv are correlation functions 
permitted by Quantum Mechanics and Local Realism, respectively. For some state we can take 
the quantum-mechanical correlation function as F. Thus we cannot falsify Local Realism if the 
inequality 

lo'' ■■■ lo"" ^QAli^l^l^ ■■■i^N)d(j)i... dcpN 
- Jo'' ■■■ Jo'' EQMi4'l, ■■■,(f>N)ELHvi4'l, ■■■,4>N)d(l)l---d(j)N (3.105) 

is always saturated, independently of the chosen LHV model. 

We will now find the necessary and sufficient condition for the violation of the inequality. The 
quantum-mechanical correlation function reads 

Eqm= EL..,„=irH....„nti4N (3.106) 

where c]^ — cos{(f)k — (z — l)7r/2). <^ H^i ^'k \ ^ mutually orthogonal vectors of 

norm tt^/^. Thus the left-hand side of (l3105l) is {Eqm, Eqm) = Ea,...n=i "^a..™ 
In Local Realism it is sufficient to consider functions in the form 

AT 

ELHv{'Pl,■■■<pN) = Y[I^'K't^^)■ (3.107) 

1=1 

The right-hand side is bounded by a specific number dependent on the correlation tensor, namely. 



{Eqm, Elhv - So'' dcpi... j^'^ dcpN Elhv Yji^,..., it, =1,2^^1 --iN Yl^ 



.Ik 



1^2 J-n-.-iN llfc=i Cfc 

< i^Tma,, (3.108) 

where T^ax is the maximal possible value of Eqm in the equatorial planes, 

Tmax ^ max EQM{(f>o, ■■■,(f>N) (3.109) 

01 ,...,0iV 

Let us now argue for this bound. Notice that in p.lOSp we deal with a sum of products of the 
integrals J^^ dcpil^^^ coscpi and J^^ dcpil^^^ sin^; with coefficients rii...i„. Thus only a projection of 
onto the subspace of normalized vectors -^cosi^i and -^sin0i is relevant. Any normaHzed 
function in that subspace can be written as 

cos0o.i-^cos</)i + sin(/)o,i^ sin0i = -^cos{(f)t ~ 0o,i). (3.110) 



As |/W(0OI = 1, one has 



/ ^ cos(</)o,^ - (pi) < (3.111) 

Jo 
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Now, since cos^i and sin are orthogonal we have 

/ d0,/Wcos0, = cos(/j,||/WII|| (3.112) 
Jo 



and 



/ d(/)i/Wsm0i =sin^,;||/WII||, (3.113) 
Jo 

where is some angle. We now can express {Eqm,Elhv) as 



iEQM,E,Hv)=n^/^ll\\^^n] E (3-114 




where = cos(iy9j — {ij — 1)tt/2). The last expression can be now bounded with l|3.109p and 
|(3JTT11 : 

{Eqm, Elhv) < ^^T^ax- (3.115) 
Thus the necessary and sufficient condition for a violation of (|3.105p is 



> 1. (3.116) 



4/ T 

It is now easy to show that the falsification of Local Realism using p.lOSp by noisy GHZ states 
l(3?7l) . p = V\GHZn){GHZn\ + (1 - V) ^^1^,^" (0 < ^ < 1), is more robust against the admixture 
of the white noise than for WWWZB inequalities. As we have explained in this Chapter, the latter 
are satisfied if, and in case of noisy GHZ states only if, J^a n=i 2 ^a...n — 1- xy-planes a noisy 
GHZ state has 2^~^ non-zero entries of the correlation tensor, all equal to ztV. The sum is equal 
to V^2^-^. From ll3ll8l) and HKl^ we obtain that for 2(2/7r)^ <V< 2-(^-i)/2 WWWZB 
inequalities are satisfied, but (|3.105p is violated. Thus the region of V, in which Local Realism is 
still possible, shrinks exponentially with N. This decay is faster for (|3.105p than in case of MAKE, 
WWWZB or WZLPZB inequalities. 

Such a derivation can generalized to inequalities, which utilize all possible apparata settings, 
when we allow di ^ 7r/2. In such a case, if all LHV-based models satisfy 



< 



sinOidOi J^^ d(j)i... sin 9^(10 N d4>N 

^E'^,^j{ei,(i)i,...,0N,4>N) 

sin 9 1 d9 1 J^^ d(f> i... sin 9i\jd9N Jq^ d<j)j\! 

'xEQM{9l,(l)l,...,9N,4>N)ELHv{6l,<t)l,---,dN,4>N) , (3.117) 



the local realistic description cannot be falsified. Similarly to the previous case of the mono- 
planar inequalities (|3.105p {6i = n/2), we introduce local bases which this time are the spheri- 
cal harmonics, in particular the three, tensor products of which span the subspace of Eqm are 

cos6'i, y^^sin^i cos0i, y^^sin^i sin0i. By an argument similar to the one described above 
they lead to a necessary and sufficient condition for p.ll7p to be violated, similar to l|3.105p : 



^ \ Z^a,.. .,n=l c 



a...n 



> 1. (3.118) 



This type of inequalities was first discussed by Zukowski |101j and Kaszlikowski and Zukowski 
|56) in papers entitled Bell Theorem Involving all Settings of Measuring Apparatus and Bell In- 
equalities Involving All Possible Measurements, respectively. Indeed, if one strictly follows p.lOSp 
or l|3.117p . observers must perform infinitively many experimental runs, at least having their ap- 
parata set in xy-planes. However, a correlation function can be reconstructed with 2^ (or 3^) 
measurement settings. Such data are enough to for p.ll6p (or p.llSp ). 
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An important thing to be noticed is that additional constraints on Local ReaHsm lead to a 
stronger version of the Bell theorem. In [64] we demand that the correlation function is invariant 
under all rotations of the Cartesian frame. This assumption leads to Bell inequalities with an 
exponential advantage over WWWZB and WZLPZB inequalities. This advantage is expressed in 
terms of the highest maximal violation ratio. One may thus ask what other constrains are helpful 
in falsifying Local Realism. 

A direct apphcation of such inequalities can be more feasible experimentally if the scalar product 
l|3.103p is defined not as an integral, but on several measurement settings per observer. Some 
applications of such inequalities are described in, e.g., |63) . 
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Chapter 4 

Critical Parameters in BJSS scheme 



4.1 Proposal of Tan, Walls, and Collett 



Most of proposed and performed Einstein-Podolsky-Rosen experiment were based on coherent 
correlations between two (or more) quantum systems. It is worthy to recall that the existence 
of entanglement emerges directly from the principle of superposition, which in general does not 
require two separate subsystems. It is thus natural to ask whether coherences of a state of an 
individual quantum system can be understood as a sing of entanglement in any situation. Can 
such states violate Bell inequalities in physically feasible experiments? The first difficulty is to 
ensure that the superposition is reaHzed in spatial degrees of freedom. Local superpositions, Hke 
an eUiptical polarization of a photon, obviously cannot provide space-Hke separation. After having 
this problem solved by an appropriate unitary transformation, we must deal with super-selection 
rules, which do not allow us to know a relative phase between components of the the states with 
different numbers of particles. In the most straight-forward attempt of a realization of a Bell- 
like experiment with a single particle the observers can at most detect the particle in one of the 
locaHzations, but they cannot reaHze any other observables, which do not commute with the first 
measurement. 

One of the first articles concerning this problem is due to Tan, Walls, and Collett [79j. They 
proposed a setup presented in Figure 4.1. It consists of three balanced beam splitter^, BSO, BSl, 
and BS2, as well as four detectors, which are able to determine number of photons detected within 
the exposure time. The action of BSO is described by 



balanced beam splitter transmits and reflects every photon with equal probabilities, regardless of the polar- 
ization. 




(4.1) 




coherent 
beam 



one 
photon 



coherent 
beam 



Figure 4.1: The scheme of Tan, Walls, and Collett. 
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and respectively for BSl and BS2: 



Cfe 

dk 



V2 I * 1 



ak 
h 



(4.2) 



with k = 1,2. X is an annihilation operator related to a spatial mode x. u and v are inputs of BSO, 
and h\, 62 are its outputs and also inputs of BSl and BS2. Their other inputs are ai and 02. The 
rest of the symbols denote outputs of BSl and BS2. Mentioning only the spatial modes we have 
assumed that all the signals will have the same frequency an polarization. The total transformation 
reahzed by the beam splitters can be written as 





( 






C2 




V d2 ) 
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V2 
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1 

i 

V2 
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0,2 
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J 



(4.3) 



Tan, Walls, and Collett first suggest the coherent states of light. 



(4.4) 



n=0 



(with a real and \n) being a normalized Fock state of n photons) to be sent to input channels Oi 
and 02 with the vacuum in channels u an f . In such a case the mean number of photons registered 
by any of the detectors is {(^^.Ck) = {dldk) = , and the mean product of the results yielded by 
an arbitrary pair of detectors equals to \a^. As expected, there are no correlations, since the two 
signals never interfere with each other. 

In the next step we introduce a single-photon signal to input u. It is easy to argue that the 
mean number of photons registered by a single detector will be increased by ;j , as the single photon 
has equal chances to propagate to any of the four detectors. Moreover, it interferes at BSl and 
BS2 with the coherent beams, which impHes the following coincidences of light intensities between 
pairs of detectors: 



{Icjc^) = {Idjd2) « (cIciC^C2) = {d\didld2) = ^(a^ + a(l + sin(6'i - 6*2))), 



7t i it , 



{Icjd^) = {Idjc2) oc {c\cidld2) = {d\dAc2) = ^(a^ + a(l - sin(6'i - 62))). 



(4.5) 
(4.6) 



The above relations suggest that the phases of the coherent beams can play a role similar to 
apparatus settings in a two-particle Bell experiments. We can now calculate the interference 
visibility according to the Michelson formula: 



V 



max{lcjd2} - min(/ci/d2) 



(4.7) 



max{lcjd2) + min(/ci/d2) 

with maxima and minima taken over ^1 and ^2- In the case of a single photon as the u- input, we 
get 

= (4.8) 

The classical analogue of this experiment is a situation, in which the u-input is not a Fock 
state, but also a coherent state Then the coincidences are 



{IM = {hjd.) <x \{a^ + a^f3\l + sin(ei - ^2)) + \f3^) 



1 



{hjd,) = {IdJc) (X -(a^ + a'(3\l - sin(0i 



02)) + \p'). 



(4.9) 
(4.10) 
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These put into ijTTl) give 



Vci = , • (4.11) 



(?) 



The relations (|4.8p and (|4.1ip are enough to argue the non-classicahty of the first situation. 

For very low values of a the term in the coincidences is neglible and thus Vqm tends to the 

unity. On the other hand, in the second case coincidences have terms j^f^'^ in the denominator, 
which are fixed with respect to a and thus the maximum of Vd = ^ is achieved for = ■ 
Any value of the visibility above ^ thus implies that the field has no classical description. This 
is not a question, however, of a possibility of local hidden variables, which would allow to model 
probabilities of results of measurements, but rather whether can the state be reconstructed using 
only coherent states, which are here seen as classical. 

The analysis is pushed further by defining a correlation function, 

- 1 < E{eM = ((^^i-4ci)(4^^-4g2)) < ^ (4^^2) 

which can be used in (|2.15p . After taking into account the coherent states of modes ai and 02, as 
well as the vacuum in v we get 

where {u^) = |(u^)|e*^. When we assume that there is one photon in u the last terms in the 
numerator and the denominator vanish, and the correlation function reads £(61,62) = '^'"j 4^ "2^" ^ ■ 

The CHSH inequality can be violated if a < Vv2 — 1. If, however, the signal sent to the channel 
u is a coherent state |/3e'^/^), we have (ii'^u) — \{u^)\ = (3'^ and {ii^'^v?) = (3^. I|4.13p function 
can be then factorized to £(61,62) = Sa^/?^ cos(6'i - ^/2) sin(e'2 + C/2)/(a^ + 2/3^)2 and hence no 
violation of a Bell inequality is possible. 



4.2 Scheme of Bjork, Jonsson, and Sanchez-Soto 



The idea of Tan, Walls, and Collett started a discussion on the non-classicality of a single 
photon |46 [ [28 t [66 t 138) . An important variant of the original proposal is due to Bjork, Jonsson and 
Sanchez-Soto (BJSS) [12]. In their scheme one is interested in probabilities of certain events, rather 
than mean values. This allows to use the CH inequality (|2.25p . which turns out to be inequivalent 
to the CHSH inequality ((2A5l) . 

Let us start with a description of the experimental setup. In its center we have a polarizing beam 
splitter (PBS), the inputs of which are denoted by a, b and outputs by c, d. It refiects horizontally 
polarized (H) light and transmits vertically polarized (V). The single photon, sent to input a, has 
a —45° polarization. Its state is thus (a^ — ay)\ft)/V^ in front of PBS and (c ^H-di)\n)/V2. The 
photon is superposed in two modes, which are distinguishable spatially, as well as with respect to 
the polarization. The other input is used to inject the coherent beam from the local oscillator with 
a +45° polarization, and the mean photon number 20^. Behind PBS the state of all four modes 
can be written as 

1^) = (e*'^^''|0,ae*'^^%l,ae''^^'^) -e*'^^<=|l,ae'"^%0,ae*'^^'')) /V2 (4.14) 

with mode ordering CH,cv,dv,dH- In the single photon modes we have used the Fock formalism, 
whereas the coherent state formalism is appHed to the other two modes. The complex phases 
appeared due to a propagation between PBS and the measuring devices for times Tc and t^, but 
they can be neglected without any loss of generality. 
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Figure 4.2: The scheme of Bjork, Jonsson, and Sanchez-Soto. 



Now, the observers shall try to detect their modes in one of the states 

1 / Ink 



-,nk,(Pk) = 



r|0,nfe) + e'^'=|l,7ife-l) 



(4.15) 



It is worthy noting that the local vacuum in both modes, |0, 0), is not covered by these states. The 
probability of having such a state at one side is e~" /2. This is one of main reasons to have a 
large, in order to minimize the occurrence of 0-photon states. 

We shall begin with demonstrating that in the perfectly realized experiment we can violate 
i|2.25p . A local projection onto | + , nt, ^/c) (+, rtfc, occurs with a probability 



P+{nk,<l)k) 



1 



(rifc-l) 



(4.16) 



1 + ^ K- ~ 1)! 

These probabilities sum up to P+{<j)k) = P+i^k, 4>k), and an interferometric relation 

P++(0c,0d) = P+(</'c)P+(0d)(l"Cos((/.,-0d)) (4.17) 

holds. 

We need to argue in what way these probabilities can be added. P+{4>k) contains a function of 
a random variable, Uk — 1, and a Poisson probability describing the distribution of this variable. 
The only parameter of this distribution is a^, equal to the mean value, the variance, and all higher 
cumulants of — 1. We can easily find these quantities by treating as a thermodynamical 
parameter and taking the partition function Z = e" : 



2 _ log Z 

~ a(Q2)2 



{{uk ~ m ~ ({uk - 1)) 

{{uk - 1)3) - 3{nk - l){{nk - 1)2) + 2{nk ~ 1) 



3 _ 3^ log Z 



(4.18) 
(4.19) 
(4.20) 



Let us first consider the variance. If is much larger than the unity, it can be neglected against 
appearing in the two terms on the left side of (|4.19p . (n^) and (rik)^ (hereafter, we take Uk ~ Uk — l 
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Figure 4.3: Accuracies of approximations l(T2T|) . l(T22|) . l(T23l) . (|i?2i|) for x = 1 (left) and a; = 0.2 



(right). * stands for (1 + a;) X]n=i l 
*** reads (1 + x) l^i 



(n-1)! ' 



** represents Z]n=i 71 



a^(n-l)! ' 



and 



appearing in averages). It is thus reasonable to assume {n^) « (fifc)^. Subsequently, we can draw a 
similar argument for the third cumulant and show that {n\) « (rik)^, etc.. Thus for any sufficiently 
smooth function we can assume {fijik)) ~ f{{nk))- In particular, let us introduce approximations 
useful in further consideration: 



^1 + % (rik-iy 



Y 1 



oo 

E 

life 



+ ^ a2(nfe-l)! 1 + x 



{a^xY^' 



oo _^ 

y i 



(4.21) 

(4.22) 
(4.23) 
(4.24) 



The role of x will be clarified below. Strictly speaking, in JMU, ([422]), l(4!23| . and ()424l) 
we demand that a^x, rather than a^, is much larger than 1. We are interested in the range 
< a: < 1. For the perfect case of a; = 1 the local probabilities obviously add up to ^, and due 
to the relation l|4.17p the inequality (|2.25p can be maximally violated, if the chosen phases are 
(j)c = 0, c^^ = 7r/2, (j)d = it/^, (j)'d = — 7r/4. The expression equals ^(\/2 — 1) > 0. 

4.3 Some Possible Imperfections [92] 



Let us now reproduce the results from [92]. Having shown that the experiment might be consid- 
ered a valid test against Local Realism, we can introduce some imperfections to our considerations. 
The first one is the coherence loss. However, since the coherent state is a superposition of infini- 
tively many Fock states, for simpHcity we make a conjecture that only the single photon part is 
affected: 

^(|0, a, 1, a) — |1, a, 0, Q;))((0, a, 1, aj — (1, a, 0, a|) 
— *■ ^(|0, a, 1, a) (0, a, 1, a| + |1, a, 0, a) (1, a, 0, aj) 

-^(|0,a,l,a)(l,a,0,a| + |1, a, 0, a)(0, a, 1, a|). (4.25) 
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We call < Z < 1 the decoherence parameter. Such an action can be justified by classicality (in 
some sense) of a coherent beam, and strong non-classicality of a single photon. 

Another imperfection we may introduce is the loss of photons during the propagation. Let us 
assume that the probability that a photon will reach the measuring devices is rj. This efifect can 
be modeled by putting beam splitters with transmittivities rj in channels c and d and tracing out 
the reflected modes. Was it possible to design a device realizing projections \+,nk,4>k}{+,nk,4>k\ 
with only one photosensitive element, the transparency of the channel would be equivalent to 
the detection efficiency of this photosensor, provided that the propagation channels are perfectly 
transparent. The single photon state is stochastically transformed into the vacuum in respective 
modes, whereas the amplitude of the coherent part is decreased by a factor y^. The whole 
transition due to both imperfections reads 

i(|0,a,l,a) - |l,a,0,a))((0,Q;, - {l,a,0,a\) p{l,T]) 

+ i(|l,a\A7,0,a^)(l,Q;^,0a^| + \0,a^,l,a^){Q,a^,l,a^\) 
- ■^(|l,aV^,0,Q;^)(0,a^, l,a07| + |0, a^, 1, a^) (1, a^, 0, a^|). 

(4.26) 

The first line of the right-hand side is the term in which the single photon is lost, the second 
are diagonal elemenst of the density matrix with single photon not lost, the last are off-diagonal 
elements. 

Using formulae (|4.2ip to the vacuum part and (|4.2ip . l|4.22p to the other part, we obtain 

(4-27) 



(^) (2-77-;cos((/.,-0d)). (4.28) 



We now choose the optimal angles to obtain — cos((/)c — (l)d) — cos(0^ — (j)d) — cos(0c ^ 4>'d) + cos((/)c — 
cf)d) — 2V2. To compute the dependence between the critical decoherence parameter and the 
critical channel transparency, above which Local Realism can be falsified, we set the left-hand side 
of l|2.25p equal to and put in the computed probabilities. The obtained equation, 

-VCRIT + '^VcRlA'^ + IcRItV2) - 31JCRIT _ „ . . 

(l + riCRirr ~ ' ^ ' 

apart from a trivial solution tjcrit = is satisfied for 

7 - 3 - 2TJCRIT + VCRIT on^ 

ICRIT - —7= • (4.30) 

2v2r?ci?/T 

Equation (|4.30p implies that if coherence is perfectly preserved, the required channel trans- 
parency is rjcRiT = 1 + ^ — 2'^/^ ~ 73.4%. This about 0.096 less than the required efficiency in 
two-photon Bell-type experiments with the maximally entangled state. 

The next step of the analysis of the BJSS scheme presented in (9^ is using the CHSH inequality 
l|2.15p . For this purpose we need to define a correlation function, which, depending on local phases, 
would take values between —1 and -1-1. We can naively do it by associating the states nfc, (/ifc) 
proposed by Bjork, Jonsson, and Sanchez-Soto with outcomes of local observables, whereas 

the states \ — ,nk,4'k) = jip^ {\^ink) — e*'^'' y^^ll, — 1)) will correspond to outcomes " — 1". 

Then the correlation function is taken as the mean value of the product of local outcomes, and 
after applying approximations (|4.2ip . I|4.22p .f l423l) . and (|4.24p we get 



E{<pc,(pd) = -V , ■ (4-31) 



3 - + 277^ -f Ar]l cos(0c - 4>d) 

After putting l|4.3ip into l|2.15p and optimizing over the angles we obtain the necessary condition 
for the violation, 

I > (4.32) 
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which impHes that the inequahty can be violated only for tj > 63.2%, significantly less than the 
value obtained from l|4.30p . Both results are correct according to the approximation. However, 
together they suggest that the analysis is not complete. The CHSH inequality is based on the 
correlation function, whereas in the CH inequality we utilize local and global probabilities. This 
fact reveals a greater generality of the latter. Thus if there exist a region of an advantage of a 
correlation function-based expression over a certain probability-based one, there must be a set 
of CH inequalities, in which at least one is violated in the region of the violation of the CHSH 
inequality. 

The CH inequality, the violation of which turns out to be more robust against the signal 
attenuation of the propagation channels than in the case shown above, is constructed with the states 
\+,nk,(j)k) being associated with the first events at each side, which enter l|2.25p and \ — ,nk,(l)k) 
with the second events. In such a case, the CH inequality reads 

p++(0c, cPd) + P+-{4'c, <f>'d) + P-+{€, M - P—i€-^ ^'d) - P+i^c) - P+iM < 0, (4.33) 

and probabilities entering (|4.33p . that have not yet been computed, are 

P—i^cM 

^{T^)\^-^os{cj,,-M), (4.34) 



p+-{(bc,M^ P-+i^c,M 

= EZ,n,=i (+' (-, na, M^^'^Piv, l)\+, nc, (/-c)!^! h, n^, <j>d) 

= Er..„.=i(-,"c>c|[^l(+,«d>d|[''V(7?,Oh,"c>c)t^l| + ,n,,0,)M 
- 2(1+^(3 - 2r; + 7f + 2lricos{^, - 0^)) 

Altogether, after optimizing over the apparata settings, l|2.25p simplifies to 

3-77 



I > 



2V2 



(4.35) 



(4.36) 



In such a case, the critical channel transparency with the perfect coherence preservation is 3— 2\/2 w 
17.2%. 

In [92] we have thus shown that a Bell test inspired by the original idea of Tan, Walls, and Collett 
|79) provides a much lower threshold on the channel transparency than the critical value of the 
similar parameter, the detection efficiency in two-qubit Bell experiments. A possible explanation is 
that the non-classicality of the state in the first case is brought in by the single photon. Coherent 
states are seen as quantum realizations of classically allowed states of the field. They are eigenstates 
of annihilation operators, which corresponds to multiplying by the intensity of the field in classical 
or semi-classical Electrodynamics. The single photon state is purely quantum, as the notion of 
a photon does not exist in classical physics. Another argument to justify this robustness is the 
following. When we loose the photon, we admix the vacuum state in the modes referring to it. The 
density matrix reduced to the modes of the single photon, in the basis {|0,0), |0, 1), |1,0), |1, 1)}, 
experiences the following transition: 



/ 





i ' 
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/ 1 
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-hi 
h 








(4.37) 



This matrix has one negative eigenvalue ofter the partial transposition, ^(1 — 77— — riY + I'^v'^), 
for any non-zero value of the product of I and rj and thus basing on the Peres-Horodecki criterion 
[671 HI] we conclude that, except for trivial cases of 77 = and I = 0, the state is always entangled. 
We thus expect that this entanglement can be revealed by Bell inequalities in a large range of the 
parameters. 
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Figure 4.4: ijcrit versus Vcrit for the two-photon experiment with a maximally entangled state, 
and versus Icrit in the scheme of Bjork, Jonsson, and Sanchez-Soto. The Clauser-Horne inequality 
is violated in the upper right region of the plot. 

All the inequalities in the BBJS scheme are never violated if / < This greatly resembles the 
admixture of the white noise (see Section 2.5). However, in this case, rather than the maximally 
mixed state, we admix a classically correlated mixture ^ (|01)(01| + |10)(01|). These correlations 
have no importance in the bases the observers perform measurements. 

One can also consider the case, in which the observers are able to predict the opaqueness of the 
channels and set the initial intensity of the coherent beam to be 2a^/r]. Then the vacuum state 
would be still admixed with a weight 1 — rj, but the approximated sums (|4.2ip . (|4.22p . (|4.23p . and 
l|4.24p are equal to ^. Thus the probabilities we use in the inequalities are equal to 

P+{(t>k) = P-{^k) = h (4.38) 
P++(,^c,0d) = P—i^cM = 3(1 - ^'7cos(./)e - M), (4.39) 
P+-{(t,c,M = P-+i(f'c,M - lil-lvcosi(t,c-M), (4.40) 

which trivially gives a condition IcritVcrit = 1/ V^- The photon loss has now the same effect as 
the decoherence. The inequalities cannot be violated for 77 < l/-\/2. 

4.4 Experiment of Hessmo et al. 

At the end of the discussion of possible non-classicality of Fock states we will consider an 
experiment performed by Hessmo et al. [43- The setup consists of a non-polarizing beam splitter, 
which receives a single photon and a coherent beam of a small intensity as input signals. The 
polarizations of the two signals are mutually orthogonal, and the state of the field behind the 
beam splitter can be written similarly to the one in the proposal of Bjork, Jonsson, and Sanchez- 
Soto, in [17] it differs only by a sign between the components. Each measuring device is built 
of a birefringent wave-plate to control the relative phases 0c and (j)d between the single photon 
and the coherent beam polarizations. Behind every wave-plate there is a polarizing beam splitter, 
which transmits the single photon with a probability and refiects it with a probability (t and 
r are taken real), thus the transmittivity and the refiectivity for the coherent beam are and 
i^, respectively. In the experiment « cos^2°. At the transmitted outputs the observers had 
single-photon detectors, whereas the reflected outputs were left unobserved. 

The probabilities of the detections are 

Pc^Pd^ ^e""''' (1 + r^ + c^^rh^), (4.41) 
the probability that the both detectors register a photon reads 

Pcd = (^2 ^ ^2^2(^ ^ ^^^(^^ _ ^ (4 42) 
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Figure 4.5: A scheme of a measuring device in the Hessmo et al. experiment. If we add a photon 
counter in a reflected mode (in a dashed box) and replace detector D with one able to detect 
exactly one photon, we would be able to perform a BJSS experiment. 



which allows to calculate the coincidence probability as 

Pcoinc — 1 Pc Pd ^ Pcd 

= 1 - e-"'''' (1 + r2 + a^r'^t^) 



_p-2ar +2rH\l + cos{(j,,-^a))). (4.43) 



Moreover, Hessmo et al. consider also the cases when the single photon does not reach any of the 
detectors and the two detections are caused by two photons from the coherent beam. The total 
probability of the coincidence is -Pcolnc — vPcoinc + (1 — '7)(1 — )^. One can minimize the 

effect of false coincidences by choosing « and a^r^ « 1. The second-order interference 
visibility found in f47] is (66 ± 2)% and (91 ± 3)% after the background correlation correction. As 
shown by Tan, Walls and Collett, both results are enough to agrue for the non-classicality of a 
singe photon. Moreover, after the background correlation correction the second-order interference 
visibility is high enough for violation of the Bell inequalities. 

Obviously, the Hessmo et al. experiment is not equivalent to the the one proposed in [12], The 
non-classicality is demonstrated not with violation of the CH inequality, but only with high values 
of the second-order interference visibility. On the other hand, to conduct the BJSS experiment one 
needs to enrich the measuring devices with detectors, which are able to reveal exact numbers of 
photons detected in reflected modes. The other necessary modification is to replace detectors that 
detect at least one photon with such that detect exactly one. As we show below, such detectors 
can be used to realize the projections utilized in the BJSS scheme under the condition that the 
second signal apart from the coherent beam contains at most one photon. Behind the PBS the 
detected state is 

\det,nk) - -j=l==el{jlr>'-^%Q), (4.44) 

where Cfc, fk are spatial modes of the PBS outputs. Since the signal enters each polarizing beam 
splitter trough only one input channel, the polarizations in spatial modes in outputs are already 
fixed. Now, let us denote the modes of the single photon and the coherent beam at each side as 
ksp and kcB, respectively. Then the action of the PBS is 



4 \ * r\( kip 



(4.45) 



so that the detectors are sensitive to states —j==^^{tklp ~\- rkljQ){—rksp + ifc^^)"'' ^|0,0). 
Finally, we include the birefringent plate and the states detected by the devices are 

\detk,nk,^k) (X (e''^'=ifc^p-t-r^cB)(-e'^'=V^fc^^ + Vift^B)"''~^|0,0). (4.46) 
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When we plug in our assumption of having at most one photon apart from the coherent beam, we 
get 

\detk,nk,<i>k) oc -rt("'=-i)V^|0,nfe)+e^'^'=(r'= -rV'=-2)|l,nfe- 1), (4.47) 

where the notation from the last section was adopted. The properly chosen position of the PBS in 
the measuring device for each individuahy would ahow to conduct the BJSS experiment. 
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Chapter 5 

Entanglement in Bulk Systems 



5.1 Introduction 

In previous Chapters we have studied entanglement between only a few quantum subsystems. 
We have assumed that the state of the system is known and can be arbitrarily manipulated. 

However, the presented methods of detecting entanglement do not seem to appropriate for 
macroscopic systems. For ensembles of the order of 10^'^ spins one may not refer to individual 
properties of every constituent of a solid sample or correlations between all of them. Instead, one 
should base one's arguments for non-classical correlations in the system on a statistical description 
of the whole ensemble. 

This is possible when the system is in a thermal equilibrium at some temperature T. Given 
the Hamiltonian H describing the system, the state is p = exp{-~ H / kT) / Z , where k is the Boltz- 
mann's constant, taken equal to 1 for all numerical purposes in this dissertation. The collective 
thermodynamical quantities are computed by differentiating the logarithm of the partition function 
Z — Tr exp(— TJ/kT) with respect to external parameters, like temperature or the magnitude of 
the magnetic field B. For example, the internal energy is equal toll ^ — g^x/^y) , the heat capacity 

C = ^ , and the magnetic susceptibility X = ^ a ^ • ■'■^ this Chapter we will show that under 
some conditions these quantities are entanglement witnesses ^80], that is their low-temperature 
values can reveal correlations of purely quantum nature. 

The problem of the detection of entanglement in a thermal equilibrium is an appealing challenge 
of the Quantum Information Theory. First, it allows to see quantum correlations as a natural 
feature of quantum systems at low temperatures. Second, thermal entanglement present in a large 
scale in solids might be important for future Quantum Computation Technology. 

5.2 Concept of Entanglement Witness 

Before we give a definition of an entanglement witness, let us recall in details the Peres- 
Horodecki criterion [671 HH] mentioned in the Sections 3.5 and 4.3. It states that under a partial 
transposition separable states always remain physical, while entangled ones can afterwards have 
a negative eigenvalue. For non-entangled states the statement is easily proved. The transposi- 
tion changes any state p into another density operator. Thus, for any two states p'^' and p^'^\ 
pl^l ® (p'^')"^ must also be positive. To finish the proof we notice that also a statistical mixture, 
J2i PiPt^ ® (pf is also a non-negative operator of trace 1, provided that Pi > and J^i = 1- 

To show that the positivity of a density operator under a partial transposition does is not 
preserved in general, it is enough to consider an entangled state of a system of dimensionality 
2 X 2 or 2 X 3. For these systems the positivity of a partially transposed state was shown to be a 
necessary and sufficient condition for its separability [50]. Let us, for example, consider the case 



49 



of a two-qubit singlet state: 
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(5.1) 



The Ti index stands for the partial transposition with respect to the first qubit. The same effect is 
obtained by transposing with respect to the second subsystem. The partially transposed matrix has 
eigenvalues ^, ^, — However, there are known states, which, despite of their entanglement, 
have a positive partial transpose |88j. 

The partial transposition is an example of a linear, positive, but not completely positive map0. 
Such an operation, mapping a state from Hilbert space onto can be generally written as 



12]n[13] 



ki,k2,h,l2 



(5.2) 



From this relation we see that with every linear map we can associate a self-adjoint matrix 
V^k^Mduh = (fci/i|[23lVl^[23l|/c2Z2)^^^l. In a short form the action of the map on the second subsys- 
tem can be written as 

e(p[12])[13] ^^2(^[12]T.p^[23])_ (5_3) 



K we assume that dimensionalities of Ti^^^ and H'^' are equal, the reverse relation can be also 
obtained with the maximally entangled state = X^i^i N)'^' 1"^^^'^ 

dimensionality of Ti'^l . Then 



, where D2 is the 



(l[2l ®e[2'3])(|^ + )[22'](^+|[22'] 

We stress that here l'^' denotes the trivial map on 
Proof: 



(5.4) 



(iW ® e)(|^+)P2'l(V^+|[22']) = (1(2] ^e) (^E.f^^^i |^)[21|^)[2'1(^'|[2](^'|[2' 

= i^EftiK>t^i(*'|[^i®ED,=i(*^i|f^''iw^[^''iN'^2)i^'^i|Zi)i^i(^2|i^i 



(5.5) 



QED. 

The mutual relation between the linear maps and the operators is known as the Jamiolkowski 
isomorphism [54]. It was also shown that 

• if e is a completely positive map, matrix W (upper indices of W are hereafter dropped) is a 
positive operator, 

• if e is a positive, but not completely positive, map, W is an entanglement witness. 

In other words, in the second case the mean value of W in any separable state is or positive, 
but there exists an entangled state, for which the expected value of W is negative. Of course, the 
threshold, above or below which entanglement is demonstrated, can be arbitrarily redefined, as in 
case of the operator. The Horodecki Family has shown |50] that any type of entanglement, 
bound [51] , distillable [10] , or directly violating a Bell inequality [6] , can be detected with a proper 
positive map, thus by a proper entanglement witness. 



positive map transforms one non-negative operator into another. A completely positive map e is such that 
1 (gi e is positive, regardless of the dimensionality of the first Hilbert space. 



50 



5.3 Internal Energy as Entanglement Witness 



One of the first significant papers concerning entanglement in solid state models is due to Wang 
and Zanardi [85j, who have considered nearest-neighbor entanglement in a spin-^ Heisenberg ring. 
The Hamiltonian of this system is given by 

/ N \ 




a . (5.6) 

Periodic boundary conditions are guarantied by an identification A^+l = l. J is a coupHng 
constant, positive for the antiferromagnetic case (AFM), and negative for ferromagnetic (FM) 
systems. The two symmetries of the system, the invariance under an arbitrary collective rotation 
of Cartesian frames of all qubits and a cyclic permutation of them require that the state of any 
two neighboring spins, p'^^' = Tr3....Arp, is 



/ Ml \ 

U2 z* 

z Us 

\ M4 / 



(5.7) 



The first noticeable fact is that due to isotropy of the system the sample is not magnetized, i.e., 
the Bloch vector of any qubit is vanishing, {a 1^1) — (a I^') = 0. Then the entries of the density 
matrix are found as 

Ul = "4 = ^(1 + J33), 
U2 = U3 = i(l - T33), 

z=^{Tn+T22 + iTi2-iT2i), (5.8) 

where Tu'S are elements of the correlation tensor (see Section 2.4). The second observation to be 
made is that, again, because of the rotational symmetry, T12 ~ T21 — 0, thus z is real. 

Wang and Zanardi have then computed the two-qubit entanglement measure introduced by 
Wootters [9l] and known as the concurrence. It is defined in a following way: let be a 

complex conjugate of a two-qubit state pl^^l. Define a matrix R = \J p^'^'^^ a^2^ p'^'^'^^* a^2^ 0-^2^ 
with non-negative eigenvalues Ai, A2, A3, A4 in the decreasing order. The concurrence is defined as 
Cone = max{0, Ai — A2 — A3 — A4}. In this specific case Cone — ^ max{0, jTn + T22I — T33 — 1}. 
Notice that Tn = T22 — — U/3NJ. The internal energy increases with temperature and is 
equal in the limit of T ^ cxi, as the thermal state is then the maximally mixed state, and the 
Hamiltonian is traceless. Therefore we conclude that the internal energy is always negative. Hence 
we have Cone = \ max {O, ~ ^)} AFM and max {O, ^ (3^ ~ ^)} ^'^^ particular, 

the nearest-neighbor concurrence at T ^ is determined by the ground-state energy Eq per spin. 
Hulten |53] found it numerically to be Eq/NJ « —1.773 for J > and N 00, whereas for FM 
systems it is known to be —1. In the latter case it we do not observe thermal entanglement at 
any temperature, whereas for the antiferromagnetic case the criterion of Wang and Zanardi fails 
above the critical temperature Tc at which U{Tc)/NJ = — 1. This temperature would be different 
for different TV. The only case, in which an AFM ring reveals no thermal entanglement, is iV = 3. 

Wang and Zanardi have linked the internal energy with entanglement using the concurrence. 
A similar, though more direct argument has been drawn by Brukner and Vedral [iQ. They have 
considered xx and xxx Heisenberg rings described by a Hamiltonian 
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with B being the external magnetic and J3 = for an xx ring, and J3 = J for an xxx 

model. Brukner and Vedral have defined the internal energy and the magnetization per site, 
U = U/N,M3 = {M3)/N = nT^^. They argue that in both cases the quantity Ll+^M^ jg ^n 
entanglement witness. 

Let us first focus on the case of J3 = J. Then the expHcit form of the Brukner- Vedral reads 



U + BM: 



3 



J 




(5.10) 



The proof of the witnessing properties was first presented by Toth, Giihne, and Cirac ^2]. Notice 
that for product states, p = p'*', each element of the sum can be bounded by 

N L J • O- ^ ^ ')PROD 

= W PROD ■{■<^^'+'^) PROD 

where the first inequality comes from the Cauchy inequality. By convexity, the same bound holds 
for separable states. 

Thus any value of (|5.10p . the modulo of which exceeds 1, immediately impHes thermal entan- 
glement. As we have said, the lowest possible value in the thermodynamical limit was found to 
be about —1.773. In fact, for the xxx models without the magnetic field the result of Brukner 
and Vedral fully agree with the calculation of the concurrence for nearest neighbors by Wang and 
Zanardi. 

The argument can be extended to more general cases of J3 ^ J and B not equal to zero. In 
particular, Brukner and Vedral consider an infinite xx ring, the analitycal solution for which was 
presented by Katsura [57]. In order to present his results we introduce a short-hand notation 
L = and K = and a function 



f{K, L, Lu) ^ V2if 2 2K^ cos 2w - AKL cos uj + L"^. (5.12) 
Now, the internal energy and the magnetization are found to be 

U = / f{K,L,Lu)taiihf{K,L,uj)duj, (5.13) 

7^ Jo 

m^-- r ^^jj^^j^t^nhf{K,:,cu)<k.. (5.14) 

These two integrals allow to determine a region of y and ^ (for J > 0) in which thermal 
entanglement is revealed by (|5.10p . The region is presented in Figure 5.1. 

Other examples of using the internal energy as an entanglement witness can be found in e.g. 
l2S[95l[8T]. 



5.4 Detecting Multipartite Entanglement with Internal En- 
ergy 

A bound for presence of multipartite entanglement in state of a Heisenberg chain was presented 
in [44] . The Authors have notices that the Hamiltonian of an xxx open chain with an odd number 
of sites can be written as 

Hxxx = J ^ W^2i-l,2i,2i+l, (5.15) 
i=l 

^We assume that B contains all the physical constants. For experimental purposes, we would need to include a 
of product of the Bohr magnetron ijlb and the gyromagnetic factor g. Then it is also necessary to remember about 
the Planck constant being a part of the spin operator. 
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Figure 5.1: (|5.10p for an infinte xx spin— ^ ring (for J > 0) as a function of the magnetic field and 
temperature. Only in the region marked with a white ,;k" the values of l|5.10p reveal entanglement 
in the thermal state. 



with W.j^k = 0**1 • a + a 1^1 • a I*^! . 

First let us show that genuine tripartite entanglement must exist if 

Proof: The bipartite entanglement bound for H^xx can be found by summing over such bound 
for three-qubit hermitian operators W2i-i.2i,2i+i, even if the do not mutually commute. We are 
interested in mean values possible only for states with at least tripartite entanglement, we allow 
two qubits to be in an entangled state. However, W2;-i,2/,2/+i is not sensitive to correlations 
between non-nearest neighbours, hence we shall assume that the state is pure and such, that 
|^[2i-i 2/ 2/+i]^ _ |^[2i-i]^|^[2fc+i 2;+2]^_ ^j^g language of the correlation tensor this means that 

Tijk = TiooTojk for any i = 1, 2, 3 and j, fc = 0, 1, 2, 3. Furthermore, we have X]i=i ^ioo ~ 1- Such 
a form of the state implies 

|(W^i23)| — \TiooToio + 12007020 + TsooToao + Ton + T022 + Toss]. (5-17) 

In the next step we choose (Tioo, T200, Tsoo) to be parallel to (Toio, T'0207 ^030) and to be 

Schmidt-decomposed, = cosa|10) — sinajOl). The non- vanishing elements of the correlation 

tensor are T03 = — T30 = cos 2a, Tn = T22 = — sin 2a, and T33 = —1. |(Wi23)| simplifies to 

|(W^i23)| = |cos2a| + |l + 2sin2a|, (5.18) 

which optimized over a gives V5+ 1. The Hamiltonian consists of (iV— l)/2 Wijk elements. Since 
we are interested in the limit of large A'', the bound obtained above is divided by 2. 

It is easy to see, however, that the bound l|5.16p is not optimal. The argument is not self- 
consistent, as we first assume, say, the third qubit to be entangled with the second, then we want 
it to be in a pure state. 

To find a better threshold, the Authors assume that the state of the ring of = 4M qubits is 
bi-factorisable, \ip) = and define two 6A^-dimensional vectors 

VI = ([1], [1 : 2], [2], [0], [5], [5 : 6], [6], [0]), (5.19) 
V2 = m, [0], [3], [3 : 4], [4], [0], [7], [N ~ 1 : N]). (5.20) 

Here H = ((afl), (a^), (a^)), [z : j] = ((^[^[^1), {al^a^^\, (afal^)) and [0] = (1, 1, 1). Then the 
mean value of the Hamiltonian is given by {Hxxx) = ^J^i ' ^2- It is important to notice that for 



53 



Non-iirfear 
wilness 



Optimal 

linear 

witness 



Figure 5.2: A non-optimal linear witness defines a hyperplane far from the convex set of separable 
states and thus can detect only strong entanglement. A hyperplane of an optimal witness has 
common point with the boundary of the set, but a non-linear witness better describes the curvature 
of of the boundary. 

any two-qubit state J^i^iC^io +^(>j + — Hence the norms l|5.19p and l|5.20p are bounded by 
'^■liHxxx)] < Jy/vi ^, JV^2 ^ < JVSN. Applying the Cauchy inequality we get that if 



the state must consist of genuine three-partite entanglement. 

In a similar fashion Giihne, Toth, and Briegel argue for the bound on the energy for bi- 
factorisable states in case of the xx ring. 

5.5 Non- linear Entanglement Witnesses 

In most of considerations presented up to this point we have discussed detecting entanglement 
with linear functions of the state, mean values of some operators. However, unlike the set of 
statistics allowed by LHV, the set of separable state is not a hull spanned between a number of 
vertices. It is rather a convex object with infinitively many extreme points (factorisable states). 
Thus we expect that its curved boundary would be better approximated by a non-linear function 
of the state than by a hyperplane. In other words, an optimal non-linear witness could detect 
entanglement in more states than a linear operator, which bounds the set of separable states (see 
Figure 5.2). A non-linear criterion of entanglement can be drawn for the gedankenexperiment of 
Einstein, Podolsky, and Rosen \S^. In that situation it can be suprising that for an entangled 
state uncertainties of a difference of positions of two particles and a sum of their momenta can 
simultaneously vanish. For product states one has A^(p[^l -|-p'^') = A^(pl^l)-|-A^(p[^]) and A^(g'[^l — 
g[2]) = A2(g[il)-HA2((j[2l). Now, since A2(gH)A2(pH) > , it is obvious that the both uncertainties 
cannot be equal to at the same time. More explicitly this argument was drawn for continuous 
systems by Reid [70] and Duan et al. [30]. 

We, however, want to focus on a version of this criterion formulated for discrete systems by 
Hofmann and Takeuchi |48] . Let us consider a spin of arbitrary magnitude I and the variances of 
its three components, first in a pure state 1-0) (ft = 1): 



It is also necessary to show that a mixedness of a state p — Pi\ipi) (ipil can only increase the 



(5.21) 



A2(5i) + A2(52) + A2(53) 

> l{l + ^l. 



(5.22) 
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variance of any Hermitian operator A: 

= E.P^\{A')^~{Af.+{{A).-{A)f 



=A^{A), >0 

>E.^.A2(A)„ (5.23) 

with {A)i — {^ilAltpi) . Hence, if we have two local spin operators S^^^ and 5^^' (the two spins 
do not need to have the same magnitude, let here us assume, however, that ^'^1 = ^'^1 = /), for 
separable states 

A2(S^W +5[21) > 2? (5.24) 
holds. On the other hand, one can have the two spins in a singlet state 15 = ^), for which 

A2(5W + 5[21) = ((S^W + 5[21)2) - + = o. (5.25) 

The argument can be straight-forward generalized for more spins. This will be the main theme of 
the next Section. 



5.6 Magnetic Susceptibility as Entanglement Witness [90] 

The criterion of Hofmann and Takeuchi can be verified for a wide class of large spin lattices 
with a thermodynamic quantity, i. e., the magnetic susceptibility. The use of this function of 
the state as a witness was first demonstrated for a specific material, copper nitrate [19], and then 
shown in its full generality by Wiesniak, Vedral, and Brukner [90] • The arguments from [19] will 
be recalled in the next Section. 

The appropriate systems to study their magnetic susceptibility from the theoretical point of 
view are those, Hamiltonians of which are invariant under arbitrary collective rotations. Let the 
system be described in the zero magnetic field by Hamiltonian Hq. In presence of the external 
field B it becomes H — Hq + B- M = H^ + B- Y^f^i S"'''. Let us consider the magnetization in 
the direction of the field B, {Mg) = -f^T^^. Using the Trotter identity, 

e^+^ = lim (e^e^Y , (5.26) 

which is true even if [A, B] ^ 0, and the fact that operators can be cyclically permuted under the 
trace, TrABC — Tr BC A, we see that independently of the Hamiltonian the magnetization along 
any direction pointed by a unit vector n is (Af,-j) = Trpn • X]t=i However, from the same two 
facts we see that the magnetic susceptibility, Xfi = ^ ggi ^ , is in general equal to 

^" " ^ (i T^iMnP''Mftp^'-)dx - , (5.27) 

where -Bj^ is the magnitude of the magnetic field projection onto the unit vector n. Only in 
cases, in which the Hamiltonian, and thus the thermal state commutes with the magnetization 
operator, the magnetic susceptibility in a direction of vector n equals to Xfi = (Mft) . Since 

we need zero-field magnetic susceptibilities with respect to three orthogonal directions to apply 
the Hofmann- Takeuchi criterion, the Hamiltonian must commute with the magnetization operators 
in all orthogonal directions. The interaction must be isotropic and no magnetic field should act 
on any spin. On the other hand, we do not have any other requirements about the lattice. It is 
irrelevant what are magnitudes of spins seated in sites of the lattice, what is its shape, and what 
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Xnls=o < (5-30) 



are exact values of coupling constants. The lattice could be one-, two- or three-dimensional, with 
or without non-nearest-neighbor interactions, or inhomogeneous. In any case we argue that 

x'Ib^o = XiIb=o + X2|b=o + X3|b=o < ~T (5-28) 

cannot be explained without entanglement. The values are summed over three orthogonal direc- 
tions. When all spins in the lattice have a common magnitude Z, (|5.28p simplifies to 

x'ls=o<^, (5.29) 
and due to the rotational symmetry it is sufficient to consider only one direction: 

Eirf 

In particular, the zero-field magnetic susceptibility at T ^ must tend to be infinite if one does 
not allow the concept of entanglement. Experimental studies of various materials give evidences, 
however, that this not the case, that is also zero-field magnetic susceptibilities of ferromagnets tend 
to with temperature. This might be due to a fact that for these materials at certain temperature 
we deal with a significant symmetry breaking, either spontaneous, or due to weak magnetic fields 
in the environment, e.g. the geomagnetic field or fields originating from the measuring apparata. 
In presence of a distinguished direction, our criterion cannot be applied. 

Let us assume that ViZ'*! — I. In [90] we have also argued that the magnetization and its 
variance, which in specific situation can be associated with the magnetic susceptibility, together 
satisfy a certain complementarity relation: 

Non-local properties Local properties 

Proof: we start with showing that 

•2, ^fNl+l\ ,,^^,2 



{M') > [^-^ j {My. (5.32) 

For this purpose we will work in the total angular momentum basis {\L,m,i)}L,m,i, where 
M^|L,m,z) = L{L + l)|L,m,i), the magnetization in the z-direction M3|L,m,z) = m\L,m,i), 
and index i is due to possible degeneracies. Since both sides of l|5.32p are rotationally invariant, we 
can choose the z-axis as the direction of the magnetization, (Mi) — (M2) — 0, (Af)^ — (AI^)^. 
Let us now define define operator K, such that K\J,m,i) — J\J,m,i). Given probabilities 
Pj.m — J2i Pj,m,i to find the system with the total angular momentum J and the magnetiza- 
tion m we notice that (M) = J2,jm Pj,m'ni < X^jm Pj.mJ =< {K). The proof is being completed 
by facts that (iW^) = {k{k + 1)) and that Nl > j implies Nl{k) > {K'^). Thus we have 

<Xn - !^(Mf > <MK + D) - ^ > ^AVO > 0^ (5^33) 

QED. 

The left-hand side of (|5.3ip has been divided into two parts. The first, as shown by Hofmann 
and Takeuchi, is strictly related to non-local properties of the state and can be positive only in 
presence of entanglement. The second part, the square of the mean magnetization depends on 
properties of individual spins. In particular, if our sample is in a singlet state, both fractions 
vanish and the inequality is saturated by the constant in the front. On the opposite, when the 
magnetization is maximal, the second part is equal to 1, and the variance of the magnetization is 
equal to Nl. The first part of the left-hand side of (|5.3ip cancels itself out and the inequality is again 
saturated. In that sense l|5.3ip describes macroscopic quantum information sharing. Information 
can be stored in either individual properties of solid constituents, or relations between them, or, 
finally, in a form, which does not saturate l|5.3ip . 

Some interesting examples of applications of (|5.28p are given in the following Section. 
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5.7 Examples 



Let us start illustrating the usefulness of the magnetic susceptibility as an macroscopic en- 
tanglement witness with repeating the argument by Brukner, Vedral, and Zeilinger [19]. It is 
drawn in a rather different fashion than in fSO], but also correctly demonstrates the principle. The 
substance under consideration is copper nitrate, Cu(N03)2- 2.5D2O, in a form of a crystal. The 
detailed structure of the crystal is given in [98] . What is the most important for this thesis is that 
the system is a dimerized spin-^ chain, that is 

HcN = E • ^ "'^'^ + • ^ • 

j 

The inner-dimer coupling constant is Ji = O.MmeV and the inter-dimer constant was found to 
be J2 = O.llmeV. The Authors assume at the beginning of the proof that, since the system is 
isotropic, the part of the magnetic susceptibility dependent on the square of the magnetization 
vanishes. Having neglected the non-nearest-neighbor contribution, they are left with the zero-field 
magnetic susceptibility in any direction bein^ 



2kT V4 

iS*" and S"^ denote an odd and an even spin (as in i|5.34p within the same dimer, respectively. Since 
both spins have magnitudes ^h, the scalar product between them is not larger in modulo than 
\fi^, provided that they are not entangled. Thus if the zero-field magnetic susceptibility implies 
presence of entanglement if 

The reasoning of Brukner, Vedral, and Zeilinger is obviously different from presented in [90j . 
On one hand, they have neglected long-distance correlations, whereas we treat contributions from 
all pairs equally. On the other, they have a priori used isotropy of the system, what has allowed 
to lower the product state bounds. It would not be suitable for our argument, which is no more 
than a thermodynamical implementation of the criterion of Hofmann and Takeuchi. The latter 
is in line just a mathematical criterion, which shall be valid for all states, not only rotationally 
invariant. The other argument for not taking this assumption initially in [90] is that we then loose 
the extensibility of the magnetic susceptibility. The extensibility is here understood as the fact 
that for a collection of not interacting, identical systems (for example, monocrystals or molecules) 
the quantity shall be proportional to their number. This is a key feature for studying bulk objects, 
as we are unable to precisely determine the population of individual systems. 

Interestingly, the bound obtained by Brukner, Vedral and Zeilinger is the same as in l|5.30p . 
despite of the described differences in the derivations. This result was subsequently confronted 
with the experimental data published in [TT] . 

Figure 5.3 presents the measured molar magnetic susceptibility of cupric nitrate measured at 
various temperatures and with crystal powdering. The solid line is an interpolation of the measured 
susceptibility and the dashed curve represents l|5.36p . The lines cross at about 5K, below which 
the presence entanglement is manifested. It is particularly interesting that the experiment was 
performed in 1963, much before the discussion on quantum information became intensive. 

Now, let us pass to theoretical applications of the criterion. We begin with analyzing results 
obtained by Xiang in |97j. Therein he computes the thermodynamical quantities of long, efficiently 
infinite, spin— ^ and spin— 1 (as well as spin-|) chains described by the Hamiltonian 

00 00 



^The aim of this paragraph is to confront the entanglement witness with the experimental data published in [11) . 
Thus, following ^9] we for the moment include all physical constants. Brukner, Vedral, and Zeilinger themselves 
have not included h in spin magnitudes. 
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Figure 5.3: The temperature dependence of magnetic susceptibility of powder cupric nitrate (tri- 
angles) and a single-crystal cupric nitrate measured at low- field parallel (open squares) and per- 
pendicular (open circles, crosses, filled circles) to the monoclinic b axis. The data and the figure 
are from Ref. [TT]. The solid curve is the theoretical curve for a dimer rescaled for the amount of 
noise estimated from the experiment. This noise is computed as the ratio of the maximal exper- 
imental value (averaged over crystal data) and the maximal theoretical value. The dashed curve 
represents the macroscopic entanglement witness (|5.36p and is rescaled in exactly the same way. 
The intersection point of this curve and the experimental one de- fines the temperature range 
(left from the intersection point) with entanglement in the substance. The critical temperature 
is around rj=^P « 5K. Note that the entanglement witness will cut the experimental curve (and 
hence yield a critical temperature) independently of a particular rescaling procedure. (The figure 
and the caption taken from |l9j1. 
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Figure 5.4: The zero-field magnetic susceptibility of an infinite Heisenberg chain of spins— ^ (left) 
and spins- 1 (right) computed in [ST]. The dotted line is our entanglement criterion l|5.30p 
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Figure 5.5: The internal energy of an infinite Heisenberg chain of spins— ^ (left) and spins- 1 (right) 
computed in [97] ■ The dotted line is the entanglement criterion (|5.10p 



The used numerical procedure, called the transfer matrix renolmalization group, is described in 
details in ^4j. Figure 5.4 presents computed zero-field magnetic susceptibilities for spin-^ and spin- 
1 chains, respectively. Dotted lines are our entanglement criteria given by (|5.30p . The lines cross 
in 5.4 at T « 1.4 and at T — 2. These results are to be compared with the critical temperatures 
found using the internal energy argument, l|5.1ip . The internal energies for spin-i and spin-1 
chains [97] are presented in Figure 5.5. If one uses the Brukner-Vedral witness rather than the 
magnetic susceptibility, one finds that entanglement is certainly present below T « 0.8 for spins— ^ 
and T « 1.4 in the case of spins- 1. The advantage of l|5.30p over (|5.1ip shall be explained with the 
fact that, unlike (|5.1ip . where only the nearest-neighbor pairs are taken into account, in l|5.3Qp we 
also include the correlations between not neighboring sites. 

Let us study lattices of spins-1 in more details. For spins— ^ the interaction, which would be 
biquadratic, rather than bilinear against spin operators, would be trivial. However, for spins of 
higher magnitudes it is an interesting isotropic interaction. We have numerically computed the 
zero-field magnetic susceptibility for the ring of six spins-1 with adjustable bilinear and biquadratic 
interaction constants. This system is described by the following Hamiltonian: 

6 ^2 

H(a) = cosa^S'W-5['+il-|-sina^(s'W-5['+il) , (5.38) 

i=l 4=1 

with the periodicity condition [7] = [1]. The magnetic susceptibility per spin is presented in Figures 
5.6 and 5.7. From the latter we see that entanglement is implied by l|5.29p between a = —jir and 
a = ^TT, and only below Tc ~ 2.66. Quantum correlations are more temperature-persistent for 
-|7r < a < than for < a < ^tt. This might be due to a fact that the biquadratic term with 
a negative factor favors states with small total angular momenta over those with large angular 

momenta. For example, the two-spin-1 operator ■ 51^1^ has the singlet state as an eigenstates 
with the corresponding eigenvalue 36, while the states of the total angular momentum 6(=2(2+l)) 
have the eigenvalue 1. In the region of the entangled ground state and positive a, the bilinear 
interaction is dominant over the biquadratic term. For a close to — fvr, on the other hand, the 
latter overcomes the former. 

Another model we would like to consider is an example of a 3-dimensional lattice. Let us 
consider the system of nine spins— ^ , eight of which are placed in vertices of a cube and interact 
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Figure 5.6: The zero-field magnetic susceptibility (summed over three orthogonal directions and 
multiplied by temperature) per spin of a Heisenberg ring of 6 spins- 1 as a function of temperature 
and a. 




Figure 5.7: The zero-field magnetic susceptibility (summed over three orthogonal directions and 
multiplied by temperature) per spin of a Heisenberg ring of 6 spins- 1 as a function of temperature 
and a plotted up to the entanglement threshold. 
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Figure 5.8: 8 spins-i forming a cube and surrounding another spin- 




Figure 5.9: The zero-field magnetic susceptibiHty per spin (summed over three orthogonal direc- 
tions and multiplied by temperature) of a cube of spins- ^ with another spin-^ in its center and 
plotted up to the entanglement threshold, 0.5 

with one another only along edges. The ninth is placed in the center of the cube and interacts 
with all others equally strongly. With the free parameter a we will this time vary the proportion 
between the constant of the coupling along edges and with the central qubit: 

H= cos a • S'Pl + S'PI • S^I^l + S^^l • S^W + S^W • (5.39) 

+ ^[5] . 5[6] + . §[7] + §[7] . 5?[8] + 5?[8] . ^[5] (5^40) 

+ 5W • ^[5] + ^[2] . ^[6] + ^[3] . 5;[7] ^ s[A] . ^[8]^ (5 

+ sina(5[9l -^^^i^I'l) . (5.42) 

In Figure 5.9 we present the computed susceptibility per spin multiplied by temperature up 
to the entanglement threshold 0.5. It can be seen that the thermal state is entangled for —0.4 < 
a < 0.43 at sufficiently low temperatures. However, the quantity presented in Figure 5.9 does not 
vanish at T = 0, hence the magnetic susceptibility would, as it seems, diverge with T — *■ 0. This 
is for the reason that given an odd number of odd-half-integer spins it is not possible to construct 
a singlet state. The ground state shall possess the rotation symmetry of the Hamiltonian at the 
expense of being mixed. However, in realistic situations this symmetry shall be broken and the 
probe shall then gain the magnetization. 

This interesting result motivated us to find a model, which in various phases has all three ways 
of the low-temperature behavior of the zero-field susceptibility at T ^ 0: 

• x'T tends to the constant above the entanglement threshold, 

• x'T tends to the non-zero constant below the threshold, 

• x'T tends to zero. 



61 




Figure 5.10: The zero-field magnetic susceptibility per spin (summed over three orthogonal di- 
rections and multiplied by temperature) of a cube of spins- 1 with another spin-i in its center 
multiplied by temperature. 

- sin a 
-cos a 
• spin-i 
o spin-1 

Figure 5.11: 6 (left) and 4 (right) spins-^ surrounding a spin-1. 

The system which satisfies this requirement is a simple modification of the previous one. Instead 
of a spin— i, we shall place a spin-1 in the middle. To save the computational resources, we shall 
also decrease the number of surrounding spins and arrange them into an octahedron (o) and a 
tetrahedron (t). The systems are drawn in Figure 5.11, and the corresponding Hamiltonians shall 
be given by 

/f(„)= cosa((5W-h5[6l).(EL2^1'l) 

+ sina^I^l -^^^i^I'] (5.43) 

for the octahedron, where = 5 for i between 1 and 6 and /t^l = 1, and 

H^t)= cosa((E-=i^''0'~^)/^ 

+ sin aS^^^ -YlUiS^^ (5-44) 

for the tetrahedron, with Zl^l = 1 and ^ being the magnitude of all other spins. These systems 
have at least 8 and 6, respectively, distinct phases with respect to a. Transitions between these 
phases can be seen in Figure 5.12 as rapid changes of values of low-temperature x' ■ In their nature 
lies an interchange of the ground state and one of the excited states. 

The susceptibilities multiplied by temperature are presented in Figures 5.13 and 5.14 (up to 
entanglement thresholds, 4 for (o) and 3 for (t), respectively). Both systems have interesting 
features. The ground state of (t) is recognized as entangled whenever the interaction between 
the spin-1 and surrounding spins is antiferromagnetic, and for (o) there exists an interval of a in 
which the thermal is satisfies our separability criterion at T close to 0, but then the magnetic 
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Figure 5.12: x' at T = 0.01 for an octahedron (left) and a tetrahedron (right) of spins-^ with a 
spin-1 in the center. Every rapid change of the value corresponds to a quantum phase transition 
of the first kind, a change of ground system. 




Figure 5.13: x'T plotted for an octahedron (left) and a tetrahedron (right) of spins-^ with a spin-1 
in the center. 



susceptibility goes under the threshold as temperature grows and low excited states significantly 
contribute to the mixture. 



5.8 Heat Capacity as Entanglement Witness |91| 



In Section 5.3 we have argued that the presence of thermal entanglement can be deduced from 
the internal energy. However, it is physically difficult, or even impossible to precisely measure 
the internal energy of an object, as it is always given up to an irrelevant additive constant. One 
can gauge the internal energy, for example, by bringing the sample close to the absolute zero 
or the infinite temperature, and then by a careful energetic balance of bringing it to a desired 
temperature. As we have already mentioned, for spin lattices and other systems described by 
traceless Hamiltonians the energy at T ^ oo is 0, whereas the ground-state energy can be calculated 
from H. Neither of the limits is attainable physically, but both can be reached sufficiently close 
to make arbitrarily good approximations. To use the zero-temperature energetic balance we need 
to know the precise form of the Hamiltonian. While heating the sample and subsequent cooling 
it back to the temperature of our interest, our sample can overcome structural phase transitions, 
Hke melting or evaporating. 

From the experimental point of view, it is much easier determine the heat capacity, the change 
of internal energy per an infinitely small change of temperature, C — ^ — ■j^A'^{H) (at a 
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Figure 5.14: x'T plotted for an octahedron (left) and a tetrahedron (right) of spins-fracl2 with a 
spin-1 in the center, plotted to entanglement thresholds, 4 and 3, respectively. 

constant volume). In this section we will show that low-temperature values of C can witness 
entanglement. In particular, we will argue that in cases, in which the ground state is entangled, 
the Third Law of Thermodynamics is related to non-separability of the state. The Law was first 
given by Nernst [65] as a statement that the entropy of the ground state is a constant dependent 
only on its mixedness. We are more interested in an alternative formulation, namely that the 
absolute zero temperature is unattainable in any finite physical process, as the heat capacity then 
tends to zero. The equivalence of these two statements is still a subject to discussions, for example 
by Langsberg [58]. A convincing argument was given by Guggenheim [l3], and based on fact that 
for finite systems, the entropy at a finite temperature, S{T) — jj" ^ dT' must also be finite. 
This puts a strong bound on the possible low-temperature behavior of C . The converse implication 
is discussed in e.g. |62j. 

The first way to demonstrate the heat capacity as an entanglement witness is direct. We 
will show that the variance of some Hamiltonian cannot vanish if one allows only product states. 
By l|5.23p it must also be true for separable states. Hence either there must appear thermal 
entanglement, or the heat capacity diverges with T — > 0. 

An example of a model with non-zero energy variance for separable states is an Ising ring in a 
transverse magnetic field, the Hamiltonian of which reads 

N N 

=±^414+^1 (5.45) 

1=1 2=1 

with the periodicity -f 1 = 1. 

Here follows the proof that no factorisable state is an ei^nstate of the Hamiltonian for B ^ 0. 
Assume, that this is not the case and the state 1-0) ~ ®i=i(a2|0) + bi\l)) is an element of the 
eigenbasis of Hjsing, with some energy E. We first need to notice that this is necessary that 
ViOi, bi ^ because of the action of the magnetic field. Then one should have 

^ {00..0\His.,ngH) ^ (10..0|g,,,„g|V>) 

(oo..o|^A) (10..0IV') ' ^ ^ 

which written explicitly is 



nf=i«2 ^uli' 



(5.47) 
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Figure 5.15: The minimal variance per site IS? {Rising) of the Hamiltonian of a transverse Ising 
ring over 2-translation invariant product states {\4>i)\4>2))'^^^'^ versus the magnetic field B {J = 1). 



with E = J2iLi The last equation simplifies to 

.bC^X -4^-^ + B = 0. (5.48) 
\aij ai 

The same relation can be now obtained for the inverse of the fraction, However, one can easily 
verify that the product of the two solutions of l|5.48p is —1. QED. 

For the sake of the numerical minimization of A'^ {Rising) we consider the the states that have 
a period of two sites: 

|0i,02> = ((cos0l/2|O)+si^0l/2|l))(cos^?2/2|O)+sin^^2/2|l)))^'^/^ (5.49) 

Such a form of the state can be justified in a following way: when we minimize the energy of 
the Ising model over factorisable states, we want to allow the local states to have the same spin 
projection onto the direction of the field and in the same time every two nearest neighbors have 
should have the opposite components in the interaction direction. This can be done with states 
l|5.49p . {a2^) is not important for any i, thus we assume that the whole vector hes in the xz- 
plane. However, calculating A"^ {Rising) we also deal with means like (cg'o-g^^') (since the state is 
factorisable, all other cross terms, like (ca'cf """^'cTg """^Vg ''"*^) will vanish). It would seem that one 
needs to consider states of the periodicity of four sites. Nevertheless, numerical calculations show 
that increasing the period does not allow to mimimize A^ {Rising) more than by using the period 
of two sites. 

Let us introduce a short-hand notation Xi = (cr[*') = sinOi and Zi — (cTg') — cos 9i. Under the 
above assumptions, the variance of the energy per spin can be written as 

A^{RIs^ng)^ TV (l + + - Sz^zf) 

± 2BN{ziZ2{xi+X2) 

+ B^^{2-xl-xl) , (5.50) 

The expression has been subjected to a numerical minimization over 9i and 02, the results of which 
are plotted in Figure 5.15. 

We can see that the variance vanishes for any case, except for B — 0. When the field-interaction 
energy becomes significantly larger than the inter-spin interaction, the ground state becomes com- 
pletely determined by the field, thus a factorisable state. Even then the variance is non-zero as 
we vary the magnetic field, not the interaction strength. Note that the same results apply to both 
FM and AFM cases of the Ising model. 

These results can be confronted with the specific heat of the system found in [57]. Therein, the 
heat capacity per spin was found to be 

CIs^ng 1 f , Uj) , 



where J = fc = 1 and f{B, uj) = ^1 - 2B cosuj + B^ . 
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1 2 3 4 5 6 
T 

Figure 5.16: The heat capacity per spin (k = 1) of a transverse antiferromagnetic Ising ring versus 
the temperature T for B = 2. The dotted curve, 0.4197/r^, is the heat capacity entanglement 
criterion. 



Figure 5.16 presents the heat capacity per spin l|5.5ip for a given magnetic field B = 2. The 
curve represents a/T"^, where a = 0.4197 is the minimal Hamiltonian variance over factorisable for 
this value of the field. Below this line the thermal state is definitely entangled. 

Thus for the Ising model, and for many other models without any rotational invariance it 
sufHcies to show that A^(_ff) does not vanish for product states to use the heat capacity as an 
entanglement witness. However, more interesting from the physical point of view models, Hke xx 
or XXX systems, fail to be applied to the same criterion. At least one of the eigenstates of their 
Hamiltonians, though not necessarily the ground state is a product state, e.g., |l...l). That is, the 
energy in the state is exactly determined, hence the variance can be set within the set of product 
states. 

In (9^ we argue basing on fundamental concepts of Thermodynamics that the low-temperature 
heat capacity, which is an easily measurable quantity, can still reveal thermal entanglement given 
the certainty that the ground state is entangled. One also needs to know the structure of the 
spectrum, if it is gap-less, or the ground-state energy is separated from the first excited state 
energy by an energy gap A. The last information to be known is the bound on the internal energy, 
below which the state is certainly entangled [HIIHIIH^, Eb, as well as the ground state energy 
C/o- 

Gap-less systems: let us first demonstrate our argument for system with a continuous bottom 
part of the spectrum per spin. Examples of such systems are half-odd integer spin chains, taken in 
the thermodynamical limit. Their internal energy per spin can be written asU = Uq + {kT)'^{cq + 
CikT + ...), where 7 > and co,ci, ... are some material-dependent constants. At low temperatures 
we may drop higher-order terms and approximate U ~ Uq + co{kT)'' . Consequently, at these 
temperatures the heat capacity shall be given as 

C = ^^^coninTy~'=^^^. (5.52) 

We can now use the bound on the internal energy to derive a similar condition for the heat capacity. 
Namely, 

C < (5.53) 

impHes entanglement within the range of the approximation. Only if the ground state is separable, 
that is C/q — Eb, (|5.53p prevents (7 ^ 00 as T approaches the absolute zero. In general, however, 
inequality l|5.53p shows that the heat capacity does not diverge at low temperatures only due to 
entanglement. The result is thus similar to the one from |90j . In both cases, the zero-temperature 
finiteness of respective thermodynamical quantities can only be explained with quantum correla- 
tions of the thermal state. 

Note that this criterion is applicable for the whole class of gap- less models in 1+1 dimension^ 
For these models the conformal field theory OSS] predicts C = (ttckT) / {3hv) , with c being the 

■^Space+Time dimensions. 
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central charge in the corresponding Virasoro algebra and v - the spin velocity. Affleck [2j has shown 
that for a chain of spins-^ c = 6^/(2 + 2Z). 

An explicit example of a gap-less model is an infinite xxx antiferromagnetic spin-i chain with 
the Hamiltonian l(537l) with 5 = 0, for which Uq = -0.443 [53J and Eb = -0.25 § Hence by 
l|5.53p the condition for quantum correlations is C < 0.386/r. On the other hand, from [97] we 
estimate that the Hnear behavior of the heat capacity (7 « |r hold for T < 0.1. In this region we 
can conclude non-classicality of the thermal state from C. 

Gapped systems: On contrary, integer antiferromgnetic spin chains, and, by definition, all 
quantum systems of a finite dimensionality belong to a universal class of gapped states. Haldane 
has argued in [45j that for these systems the low-temperature heat capacity per spin reads 

hTV ( A 



C = c'^^J exp^--J, (5.54) 

where c',i5 are again material-dependent constants. The internal energy can hence be found as 
C7(r) = C/q + \l G{T')dT' = f/o + c'|^r(-<5-l,^). Now, according to [I] we approximate 
V{a,x) ~ e^'-^x'^^^ for large x and obtain 

[7 = i/o + c'KT*+2e"^/A. (5.55) 
Finally, we again use the bound for separable states and argue that 

< — ^2 — (^-^^^ 

can be only due to entanglement, of course, within the range of the approximation. Let us here 
comment that the second approximation, kT « A, is usually much weaker than for the gap-less 
systems. 

To exemplify this version of the argument, we consider an infinite spin-1 one-dimensional xxx 
antiferromagnet with J = 1, for which c' = A^/^/V^^, A = 0.411, Uq = -1.401 [Ml p] and 
Eb = — 1 (the proof is a simple extension of l|5.1ip to spins-1). Thus all values of C below 
0.165/r^ ensure us about thermal entanglement in the ring at T « 0.4. 



^Now we use spin operators rather than Pauli matrices, hence the values are divided by 4. 
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Chapter 6 



Summary 



Entangled states [73] are not only a counter-intuitive feature of Quantum Mechanics [34], but 
also a key resource for many quantum information processing tasks, like the teleportation [9], 
Quantum Cryptographic Key Distribution [8l!35l[75], or Quantum Computation algorithms, e.g., 
[121 [771126]. Some of this applications strongly rely on the discrepancy between Quantum Mechanics 
and Local ReaHsm. 

For this reason this dissertation has been devoted to the problem of detecting entanglement 
in various physical systems. In general, this problem can be approached from two different per- 
spectives. First, we may want to confirm whether a state directly violates Local ReaHsm. In some 
models a direct application of a Bell inequality might be not feasible due to a high complexity of 
the physical system. In such a case, we can apply an entanglement witness. Quantum correlations 
confirmed by such a criterion not necessarily falsify Local Realism, but can be post-processed [10] 
in order to do so. While Chapter 2 aims to present a historical outline of the Bell theorem. Chap- 
ters 3 and 4 present more recent results, and in Chapter 5 we focus on discussing entanglement 
witnesses, particularly those, which can be deduced from thermodynamical properties of solids. 

We have begun with presenting the apparent paradox noticed by Einstein, Podolsky, and Rosen 
[34], who had considered two distant particles in a quantumly correlated state. Basing on Relativity, 
the Authors assumed that no useful information can propagate between two remote observers with 
a superluminal speed. They also had a belief of Realism, that is that results of all possible 
experiments are predefined. However, the position of one particle from the EPR pair could be 
determined by the measuring the position of the other, regardless of the distance between the 
two particles. Similarly, the momentum of one of the particles can be known by measuring the 
momentum of the other without any disturbance of the first one. Thus both the position and the 
momentum are seemingly elements of the reality. However, the quantum-mechanical Heisenberg 
uncertainty principle forbids to know precisely the values of two non-commuting observables, Hke 
the position and the momentum. This brings the alternative of two possibilities. Either Quantum 
Mechanics is unable to provide the most complete description of the reality, or these two quantities 
cannot be simultaneously elements of the same physical reality. 

To answer this question Bohr |16j has referred to his own concept of complementarity. According 
to him. Quantum Mechanics is a complete theory. Nevertheless, according to him, the act of 
measurement disturbs the system in a uncontrollable way, hence results of future measurements 
are irreversibly disturbed. The problem also lies in the distinction between a quantum system and 
a measuring apparatus. Simply, as Bohr states, a measurement of the position corresponds to a 
different physical situation that a measurement of the momentum. 

The possibility of completing Quantum Mechanics with additional parameters, which would 
determine results of measurement, yet remain unknown to observers, was considered in 1952 by 
Bohm [m dl]. He first proposed his own interpretation of Quantum Mechanics, based on an 
assumption, that while the position of a particle is stochastically described by a square of moduli 
of the wave-function, its momentum is given by a gradient of the phase function. Then stationary 
states, which have their wave-functions real, do not express any motion of the particle. In such 
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a case the particle can be seen as an entity from an ensemble of identical systems. Elements 
might have a precisely defined position and the wave-function would rather now describe the 
density within the ensemble. Thus both the position and the momentum of the particle could be 
simultaneously measured with arbitrarily high precisions. In more general cases Bohm suggests 
that the Schrodinger equation could have an inhomogeneous term, which would allow to introduce 
additional parameters and build a new measurement theory. 

However, Bell's seminal paper ^ has shown that if such hidden parameters exist, they cannot 
be of a local character. Bell showed that Quantum Mechanics, or, in particular, entanglement, is 
in a sharp disagreement to Local Realism. The discrepancy was demonstrated by inequality l|2.12p . 
the violation of which was later tested in various experiments [H [86l [71] . None of the experiments 
was fully convincing, however. In many cases, the violation of a Bell inequality has been found 
related to quantum-over-classical advantages in various quantum information processing tasks, e.g., 
cryptographic key distribution [75] or communication complexity problems [21] . 

Greenberger, Horne, and Zeilinger showed that the violation of Local Realism can be demon- 
strated, even in a stronger version than for two qubits, for larger systems. Initially formulated 
without inequalities, their argument was rephrased by Mermin [61] with help of a series of statis- 
tical expressions with local realistic bounds. Gradual increase of the arbitrarity of measurements 
performed by obsevers was presented subsequently by ArdehaH [3], and Belinskii and Klyshko |5]. 
The further generalizations were made by Werner and Wolf [88J, Weinfurter and Zukowski [87] . 
Zukowski and Brukner |103j . Wu and Zong [96], and Laskowski, Paterek, Zukowski, and Brukner 
[59] . On the other hand, there were many attempts to find Bell inequalities by finding hyperplanes 
containing faces of a certain convex hull, e.g. [36l EH [75] . A convenient way of analyzing such a hull 
for correlations between many qubits was presented recently by Zukowki [102] and developed and 
applied directly to TV = 3 by Wiesniak, Badzi^g, and Zukowski in [93]. We focused on the simplest 



non-trivial case three observables per site. The found inequalities are special cases of p.38p . We 
also show that a Bell inequality, in which all three observers choose between three observables. 

In Section 3.12 we consider effects of taking the rotational invariance of the correlation function 
as an addition constraint on theories with LEV [6^. As it appears, for GHZ states the inequalities 
which utilize this assumption are exponentially more robust against the white noise admixture than 
MAKB, WWWZB, or WZLPZB inequalities. Expressions derived in [64] utilize a continuum of 
possible local observables, nevertheless the amount of violation can by computed from the condition 
l|3.116p with the knowledge of only 2^ mean values. 

Another trend in the studies on the Bell theorem is a discussion on realizations of a Bell test 
with a single photon, rather than a set of entangled particles, as a carrier of non-classicality. This 
idea was initially proposed by Tan, Walls, and Collett [79]. One of our results [92] is showing that 
the realization of Tan-Walls-Collett-type proposal described by Bjork, Jonsson, and Sanchez-Soto 
[T2] can be suprizingly robust against photon loss, which in two-photon experiments is equivalent 
to the detection efSciency. The required detection efficiency in standard Bell experimental set- 
ups was found by Garg and Mermin [37] to be 82.8% for a maximally entangled state, whereas 
Eberhard [31J has shown that in passing from the singlet state to the limit of a factorisable state 
the efficiency can be lowered to 66.7%. We have obtained that it is necessary to preserve at least 
17.2% of photons from the source in the BJSS scheme. Interestingly, this proposal is the case of a 
great advantage of the CH inequality l|2.25p over CHSH (|2.15p . Another interesting feature of the 
schehme is that the optimal realization of the BJSS experiments requires a strong coherent beam, 
unlike the Tan-Walls-Collett scheme, which works out best in the limit of a very weak coherent 
part of the field. The considered decoherence model turned out to be analogous in its role to the 
depolarizing channel acting on a pair of entangled photons. One should mention, however, that 
with the current state-of-the-art technology we are not able perform the experiment, as it requires 
detectors, which distingush between different photon numbers. Nevertheless, there are no reasons 
to assume that the scheme will never be experimentally feasible. 

Chapter 5 focuses on the problem of detecting thermal entanglement is solids. Conducting 
a Bell-type experiment on a set of about 10^"^ spins does not seem possible. Thus one detects 
entanglement in bulk bodies with witnesses, operators or state functions, which can take some 
values only for entangled state. Entanglement detected by witnesses can falsify Local Realism, but 
can also fall into the category of bound entanglement, which cannot be even distilled to a useful 
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form [511 [To]. We begin with explaining the concept of an entanglement witness [671 HSl [80] and 
showing that for specific systems low-temperature values of the internal energy cannot be explained 
only with separable states [85 } (82 1 [TS l IST] . 

Nevertheless, we stress that it is a technical problem to learn the value of the internal energy 
at a given temperature. The first solution suggested by Wie sniak, Vedral, and Brukner [90] is to 
measure the magnetic susceptibility, whenever the Hamiltonian is rotationally invariant. This is a 
well-known routine in solid-state physics. We have proven that for these models the susceptibility 
converges to 0, or, more generally, a finite value at T ^ only because of entanglement of the 
thermal state. Notably, in contrast to the internal energy, the bound on which for separable states 
must be computed in each case individually, the magnetic response of the material is the most 
universal thermodynamical entanglement witness known up to date. As the coupling constants do 
not enter the expression, the quantity, thus its entanglement witnessing properties are independent 
of an exact form of the Hamiltonian. Our criterion can be as well applied to a simple spin-^ xxx 
chain as to dimerized, frustrated, or higher-dimensional lattices. Even lattices with spins of various 
magnitudes can be tested for entanglement with our method. The other feature of the magnetic 
susceptibility is its advantage over the internal energy in a similar role. On the other hand, however, 
it is not possible to detect genuine multipartite entanglement with x, as it was shown possible for 

U by Giihne, Toth, and Briegel [44]. It is sufficient to recall that |*) = ( ^(|01) - |10))j is 
an A^-partite state of the total angular momentum 0, which is only bipartite-entangled. 

In [90] we have also presented a non-trivial complementarity relation between the magnetization 
and its variance (which in the absence of the magnetic field can be deduced from the susceptibility) , 
which expresses macroscopic quantum information sharing. The left hand side of (|5.3ip has been 
divided into two terms and only one of them can be maximized at a time. The first is dependent 
on pairwise correlations between qubits, while other is based on their local and individual prop- 
erties. Thus the quantum information can be stored either in the magnetization, or the pairwise 
correlations, or in a form irrelevant to l|5.3ip 

Finally, we show in [91] that whenever there is a certainty that the ground state is entangled, 
also the heat capacity can serve as an entanglement witness. Without purely quantum correlations 
C would diverge to infinity at low temperatures. However, the low-temperature behavior of the heat 
capacity is under requirements of the Third Law of Thermodynamics [65] . We have also based our 
argument on the universality of thermodynajxiical quantities. In this context the universality means 
that the thermodynamical functions have similar forms for various physical systems. Thermal 
entanglement can be revealed by the heat capacity for both general universality classes, i.e., systems 
with gap-less and gapped spectra. We, however require the knowledge of the ground-state energy 
and the bound for separable states. For this reason our considerations cannot be treated as a 
derivation of new argument, but rather a version of arguments based on the internal energy. 

In conclusions, we hope that our research has brought an impact in detecting entanglement in 
various physical systems. We have presented new Bell inequalities, methods of their derivation and 
conditions for their violation [64 l (92 1 [93], The other problem addressed by Wiesniak, Vedral, and 
Brukner [90l [91] was thermal entanglement in bulk systems. In one of the papers we have derived 
thus far the most general macroscopic entanglement witness, the magnetic susceptibility. The other 
article associates quantum correlation of the ground state with the Third Law of Thermodynamics. 
We believe that all these results have enlarged our knowledge about the role of entanglement in 
Nature. We hope that in future they will also find more specific applications in processing quantum 
information. 
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